For Reference 


NOT TO BE TAKEN FROM THIS ROOM 


Gx swpnis 
UNIOMASUTATIS 
HABERTAEASIS 


Digitized by the Internet Archive 
in 2023 with funding from 
University of Alberta Library 


https://archive.org/details/Anwar1972 


% ro fe | 
' oa : ce : i? ; 
‘ 7 | i aviey "Fr. 


ale: . : yn wut : 
nin ay af ; is vo 


ORT ie . Pha Ake 


: : : | | of ina aa il 
> ae ine am 
Ores ri ah ey 


ere | ae ca . 
A NIA Ok Sad LY, ae 


THE UNIVERSITY OF ALBERTA 


EFFECTS OF NONUNIFORM MAGNETIC AND ELECTRIC FIELDS 
| IN MAGNETOHYDRODYNAMIC SLIDER BEARINGS 


BY 


(C) MOHAMMAD IQBAL ANWAR 


A THESIS 
SUBMITTED TO THE FACULTY OF GRADUATE STUDIES AND RESEARCH 
IN PARTIAL FULLFILMENT OF THE REQUIREMENTS FOR THE DEGREE 
OF DOCTOR OF PHILOSOPHY 


DEPARTMENT OF MECHANICAL ENGINEERING 
EDMONTON, ALBERTA 


Ale Gpem lie dee 


HOMAISNIA GUA ealduTe ras 3 30 yriliona aut ‘OT Garrimeue 


339099 FHT 0; 2THAMARTUDIA ANT 20 ne eens 
rHsoeoNNa 9O\OTo00 30 


lee 


UNIVERSITY OF ALBERTA 
FACULTY OF GRADUATE STUDIES AND RESEARCH 


The undersigned certify that they have read, and recommend 
to the Faculty of Graduate Studies and Research for acceptance, a 
thesis entitled "EFFECTS OF NONUNIFORM MAGNETIC AND ELECTRIC FIELDS 
IN MAGNETOHYDRODYNAMIC SLIDER BEARINGS" submitted by MOHAMMAD IQBAL 
ANWAR in partial fullfilment of the requirements for the degree of 


Doctor of Philosophy. 


ABSTRACT 


An analytical study is made of the magnetohydrodynamic effects 
in a liquid-metal-lubricated slider bearing where a magnetic field is 
applied perpendicularly to the bearing surfaces. The analyses are 
carried out for an open circuit condition, and for the case when the 
electrical power from an external source was supplied to the bearing. 

For an open circuit condition, the effects of inertia and 
the nonuniformly applied magnetic fields are investigated. The re- 
sults indicate that the effects of inertia on load capacity decreases 
with the increase of Hartmann number, M, and become almost negligible 
at high M. The highest bearing load capacity is obtained at a given 
Hartmann number when the film thickness ratio is optimized. The non- 
uniform magnetic fields i.e., increasing magnetic field, give higher 
load capacity than the comparable uniform magnetic field. The influence 
of surface conductivity on load capacity will depend upon the distri- 
bution of the magnetic field. For uniform and increasing magnetic 
fields the bearing surfaces should be insulated, while for a step-type 
distribution the stator should be a conductor. 

For the case of externally supplied electric power, a parallel 
plate slider bearing with uniformly applied magnetic field was analyzed. 
The nonuniform electric field in the bearing was created by connecting 


the segmented side electrodes to the power supply. The results suggest 
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that it would be advantageous if the electrical power is applied near 


the inlet portion of the bearing. 
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CHAPTER I 
STATEMENT OF THE PROBLEM AND 
REVIEW OF RELEVANT LITERATURE 


1.1 Introduction 

In view of the increasing use of liquid metals in engineering 
and industry, considerable attention is being paid to various magneto- 
hydrodynamic machines and devices, including MHD bearings. The design 
of bearings employing ordinary lubricants for units where the interior 
cavities are filled with liquid metal presents considerable difficulty. 
The major difficulties are due to the limitations of ordinary lubri- 
cants and the need for complicated seals. In such cases, the use of 
the working fluid itself as the lubricant, will introduce a significant 
simplification. From hydrodynamic and thermodynamic viewpoints, liquid 
metals offer several advantages over conventional lubricants. An 
important advantage is the possibility of operation at high temperatures. 
In addition, the high thermal conductivity of liquid metals enables the 
heat generated by viscous action to be readily conducted away from the 
source of generation. The resulting effect is a tendency toward uni- 
formity of temperature and viscosity within the lubricant film. The 
properties of liquid metals which can effect the performance of bearings 
adversely are low viscosity and corrosivity. 


In evaluating the possible merits of various types of bearings 
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for liquid-metal applications, it is necessary to consider the charac- 
teristics of each type of bearing in relation to the fluid properties. 
It is shown by Bisson and Anderson [1] that one important requirement 
of rolling-element bearings is an absolute compatability of the race 
and the rolling-element materials with the lubricant. This require- 
ment for absolute material compatibility with the lubricating fluid, 
makes it impractical to use the highly corrosive liquid metals as 
lubricants in rolling-element bearings. 

However, the hydrodynamic and hydrostatic type of bearings 
can tolerate some corrosion and surface damage while still maintaining 
their function. Thus, the requirements for material compatability are 
not as rigorous as those for rolling-element bearings. In hydrodynamic 
and hydrostatic bearings, a continuous film maintains separation of 
the surfaces in relative motion, and the pressure required to support 
the load is either supplied from an external source (hydrostatic 
bearing) or is generated within the bearing itself (hydrodynamic). The 
hydrostatic bearing is most suitable for supporting loads with no 
rotation or with low angular speeds. It has the advantage that for a 
fixed pressure, the load capacity is independent of the lubricant 
viscosity. Therefore, the performance of liquid-metal lubricated 
hydrostatic bearings is neither seriously affected by the high corrosivity 
nor the low viscosity of liquid metals. 

On the other hand, viscosity is the most important property 
of a lubricant in hydrodynamic bearings. The load capacity of an 


hydrodynamic bearing is dependent on the viscosity of the lubricant, 
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i.e.,the higher the viscosity of the lubricant, the larger the load 
capacity. The viscosities of the proposed liquid metals are approxi- 
mately one to two orders of magnitude less than those of ordinary 
lubricating oils at the same temperature. Therefore, the reductions 
in load capacity of the same order can be expected when these metals 
are used as lubricants. To overcome the problem of reduced load 
capacity the application of magnetic fields has been proposed. This 
has motivated the study of magnetohydrodynamics in lubrication. 

Magnetohydrodynamics is the study of the motion of an 
electrically conducting fluid in the presence of a magnetic field. 
Electric currents induced in the fluid as a result of its motion 
modify the field. At the same time, the flow in the magnetic field 
produces mechanical forces which modify the motion. The current 
engineering interest in electromagnetic fluid interaction phenomena 
has been mainly due to the concept of MHD power generation. 

Other engineering devices such as the MHD pump and meter, MHD coupler 
and bearing, etc., have been the outcome of the application of the 
theory. 

In general, the bearing configurations possible in hydro- 
dynamic lubrication may also be considered for the MHD bearings. How- 
ever, the combination of electromagnetic and hydrodynamic effects will 
further increase the number of possible arrangements available for 
study. Therefore, in the present study the analysis will be confined 


to a plane slider bearing. 
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1.2 Review of the Associated Literature 

The magnetohydrodynamic effects were first demonstrated by 
Hartmann and Lazarus [2] in their experiments in 1937. They showed 
that for a typical set of data using mercury flowing in a channel, 
the electromagnetic pressure gradient may be as high as ten times the 
hydrodynamic pressure gradient at the same flow rate. Later, Synder [3] 
in his theoretical analysis showed that the effects noted by Hartmann 
may be utilized to increase the load carrying capacity of a bearing 
operating with an electrically conducting lubricant. This created the 
interest of other authors to study the magnetic effects in the field 
of lubrication. Hence, the analyses were extended toward optimizing 
the design variables which will result in the maximum load capacity. 

Osterle and Young [4] considered a step-type slider bearing 
with a uniformly applied magnetic field. They concluded that for maxi- 
mum load capacity of the bearing, the film thickness should be a step- 
function, while the bearing surfaces should be insulators. 

Hughes [5] analyzed an inclined slider bearing under general 
electrical loading conditions. The magnetic field was applied normal 
to the bearing surfaces which were assumed to be perfect insulators. 
He found that significant increase in load capacity can be achieved 
even at low Hartmann numbers if the power is supplied to the bearing 
system from an external source. Hughes [6] also extended his analysis 
to finite width step-type slider bearing. 


Kuzma [7] considered a parallel plate slider bearing with 
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stator as a perfect conductor and slider as an insulator. He showed 
that this bearing, with nonuniform applied magnetic field, has a 
greater load capacity than a step-type bearing subject to uniform mag- 
netic field. Shukla [8] pointed out that the advantages noted by 
Kuzma are more dependent on the nonuniform conductivity of bearing 
surfaces than on the nonuniformly applied magnetic field. He showed 
that the uniform magnetic field is more advantageous provided that the 
conductivity of a bearing surface is a step-type function. His results 
for such a case indicate that for a maximum load capacity the bearing 
film thickness should be a step-function for Hartmann number, M < 5.5, 
and parallel for M > 5.5. 

Prakash [9] carried out the analysis of a composite slider 
bearing in the presence of an applied magnetic field. He concluded 
that an MHD composite bearing does not always give an increase in 
load capacity as compared to an equivalent inclined slider bearing. 

The review of literature further reveals that there has 
been very few experimental attempts to verify the theoretical results 
of MHD bearings. This may be attributed to the complications involved 
in designing the apparatus and in the difficulty of measuring the 
dependent variables. 

Maki, Kuzma and Donnelly [10] carried out the investigation 
for a step-type hydrodynamic thrust bearing theoretically and experi- 
mentally. The magnetic field was applied normal to the bearing surfaces 
and at the same time the side electrodes were connected to an external 


power supply such that there was a radial flow of the current. Using 
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mercury, the pressure and voltage distributions, and the torque at the 

slider were measured. The authors obtained good agreement between 

the theory and the experimental results. It should be pointed out, 

however, that the data were obtained for low speed and very low 

Hartmann numbers. Nevertheless, the experiment revealed that the 

increase in pressure due to magnetic and electric effects is possible. 
Other cited experiments have been performed for a MHD 


hydrostatic bearing [11], and a journal bearing [12]. 


1.3 Statement of The Problem 

The review of literature indicates that the MHD bearing 
analyses have been primarily confined to bearings where the applied 
magnetic field has been assumed uniform in part or full length of the 
bearing. Secondly, where the external power was also applied to the 
bearing to further improve the load capacity, it was assumed that the 
electric field was uniform throughout the length of the bearing. 

For the above cases where the magnetic and electric fields 
were assumed uniform, the problem was simplified from the mathematical 
viewpoint. Therefore, in most of the cases an analytical solution 
was possible. 

From hydrodynamic theory of lubrication it may be recalled 
that the pressure rise in the bearing is due to the change in the film 
thickness. This change, in practice, is brought about by the geometry 
of the bearing. Therefore a number of bearing configurations are 
possible. Typical examples of bearings in use are plane tapered, 


step-type, journal, and composite bearings. However, when a magnetic field 
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is applied to such a bearing, and if the lubricant is an electrically 
conducting fluid, then there results an electromagnetic body force, in 
addition to the hydrodynamic forces in the fluid film. The combination 
of these effects i.e., hydrodynamic and magnetic, changes the flow in 
the bearing, thereby causing modifications in the pressure distribution. 
Such effects may be utilized to increase the load capacity of the 
bearing. 

In general, the magnitude and the distribution of the body 
force produced in a particular type of bearing depends upon: 

(i) the applied magnetic field, 
(ii) the bearing surface properties, 
(iii) the external circuit conditions. 

The effects of bearing shapes and the properties of bearing 
surfaces have been investigated to some extent, but the effects of the 
applied magnetic, and applied electric fields have been confined only 
to uniformly applied fields. Therefore, it is considered important to 
investigate the effects of nonuniformly applied magnetic and electric 
Fields. 

From the literature it is noted that the inertia terms in 
the equations of motion have been neglected. However, in the case of 
a bearing operating at high speeds and utilizing a lower viscosity 
lubricant such as liquid metals, the contribution of inertia terms 
may become significant. 

In summary the present study concerns: 


(1) Nonuniform magnetic field effects in MHD slider bearings. 
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(ii) Inertia effects on load capacity of MHD bearings. 
(iii) Nonuniform electric field effects in MHD slider bearings. 

The outline of the present study stated above shall be pre- 
sented in the following format. 

The governing equations for the study of MHD lubrication are 
developed in Chapter II. These equations are then simplified for the 
particular case under consideration in the subsequent chapters. 

In Chapter III, the effects of nonuniformly applied magnetic 
fields and inertia effects are described. The analysis of nonuni- 
formly applied magnetic fields for an arbitrary bearing profile are 
presented in Chapter IV. Chapter V is devoted to studying the effects 
of step-type applied electric field. The conclusions and suggestions 


for future investigations are discussed in Chapter VI. 
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CHAPTER II 
THE DEVELOPMENT OF GOVERNING 
EQUATIONS FOR MHD LUBRICATION 


2.1 The Basic Equations 


The theory of magnetohydrodynamic lubrication is the same as 
conventional hydrodynamic lubrication, except that the interaction be- 
tween the conducting fluid and the electromagnetic fields must be in- 
cluded. As a result the number of equations to be dealt with increases. 
This interaction phenomena couples the electromagnetic and fluid vari- 
ables, thus creating additional difficulties for the solution of the 
problem. 

The basic equations for the study of magnetohydrodynamics 
are fluid equations, Maxwell's equations, Ohm's law, and constitutive 


equations. 


2.2 Simplifying Assumptions 


Prior to the formulation of the problem in MHD lubrication, 
several assumptions are made in order to simplify the basic equations. 
The assumptions to simplify the fluid mechanic equations are: 

(1) The fluid is Newtonian. 
(2) The flow is laminar. 


(3) The lubricant is incompressible. 
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(4) The lubricant viscosity is constant. 

(5) Elastic distortions are negligible. 

(6) Temperature variations in the fluid film are neglected. 
(7) Steady state has been reached. 


In addition to the fluid mechanic assumptions other assump- 
tions of electromagnetic nature are made. These assumptions include: 
(8) The conductivity of the liquid metal is constant. 
(9) There is no charge accumulation. 
(10) The liquid metal is homogeneous and non-magnetic. 
(11) Displacement currents are negligible with respect to con- 
duction currents. 
(12) All velocities are very small compared to the velocity of 
light. 
(13) Magnetostatic forces are negligible. This is justified 
since liquid metals are non-magnetic. 
Considering the assumptions above the equations of magneto- 
hydrodynamic may be written in a general form such as: 
Navier-Stokes Equations 


o[(V . v)V] = - vp + uv" 
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Vere (a) 


Continuity Equation 


vV.V=0 (2.2) 


Maxwell's field equation 
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Ohm's Law 
J = o(E + Vx B) (2.4) 
and constitutive equation 
B=u H (2.5) 


The above set of equations, in general, describe the action 
of magnetohydrodynamic effects in liquid-metal lubricated hydrodynamic 
bearings. To analyze a given lubrication problem, equations (2.1) 
through (2.4) must be solved simultaneously subject to a given set of 
fluid mechanic and electromagnetic boundary conditions. In this 
general form, however, these equations are not tractable and additional 


simplification must be made for any particular bearing configuration. 


2.3 The Governing Equations 


The bearing considered in present work is shown in Fig. 2.1. 
The coordinate system is fixed with respect to the stationary surface 
of the bearing. The magnetic field is applied perpendicularly to the 
bearing surfaces, and it is assumed to have an arbitrary distribution 
along the length of bearing. The slider is assumed to be made of in- 
sulating material while the stator may be of insulating or perfectly 
conducting material. 

For the bearing under consideration it is now possible to 
make further simplifying assumptions. These are as follows: 


(1) Lubricating film is very thin, Pe. 


(2.3) Hy =o 
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(2) Order of magnitude consideration may be used to simplify 
the equations. 
(3) The magnetic field induced by currents is assumed to be 
negligible with respect to the applied magnetic field. (Rm << 1). 
Under these assumptions equations (2.1) through (2.5) reduce 


to: 


ou ou OUs cea, OD: dU 
e(u aX ap WY ay + W aS) AN ae = us BY (276) 
a, 
On=tsr (2-7) 
(u DL ay Le pe an eh. ob aw + J.B (2.8) 
PA Ox ay az Mi at x "y 
ou aV OW _ 
Sxundoy phaztaly (2.9) 
uF = o(E, -W B.) Zen) 
J, = o(E,+ UB.) (220) 


The above magnetohydrodynamic equations have been simplified 
to a great extent, however, the equations are difficult to solve. 
In the following Chapters III through V, these equations will be 
further simplified by making additional assumptions for the particular 
model under consideration. The boundary conditions shall be specified 


and the solution will be presented. 
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CHAPTER III 
NONUNIFORM MAGNETIC FIELD EFFECTS 
IN MHD SLIDER BEARINGS” 


3.1 Introduction 


In this chapter, the effects of the nonuniformly applied 
magnetic fields are investigated. In addition, the inertia effects 
are also analyzed. In the analysis the bearing is assumed to be 


infinitely long in z-direction. 


3.2 The Governing Equations and the Associated Boundary Conditions 


For infinite bearing model, the equations (2.6) through 
(2.11) are further simplified under the following assumptions: 
(1) No fluid flow in z-direction. 
(2) The bearing and the slider surfaces are made of perfectly 
insulating material. 
(3) The side electrodes to be perfect conductors. 


Hence the governing equations may be reduced to the following 


form: 


UN eg op! onu 
7 tu ay o(E,+uB,)B, ee 


a portion of this chapter was presented in the ASLE/ASME Lubrication 
Conference, Pittsburgh, Pennsylvania, October 5-7, 1971. Later 
published in Journal of Lubrication Technology, Trans. ASME, Series 
Feu VOleeo4e, NOs, eval. 972. 
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with the following boundary conditions 
u(x,0) = U 
u(x,h) = v(x,0) = v(xsh) = 0 (374) 


p(0) = p(L) = 0 


As the electrodes are assumed to be very good conductors compared to 
the fluid, we have from, V x E = 0, that ES is a constant (as the vari- 


ation in z direction has been neglected) which will be determined 


later. 
Now since V x E = 0, a terminal potential can be defined 
as: 
te 
has | E dz (S25) 
fe) 


3.3 Equations in Dimensionless Form 


The governing equations may be written in dimensionless 


form by introducing the following dimensionless quantities. 
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Defining the stream function such that 


cect A ch 
oy ox 
Equation (3.7) becomes 
ES a oe = 
* — = aes 
pe(V ov. Wey. Pee wee Yo ug E (3.10) 
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The stream function may be expanded, after Synder [13] 
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f,(n) se (3511) 
n=0 


where n = hey/h . Oe h/h.. Substituting (3.11) into (3.10) with 
K = h,/Lh' and h' = dh/dx = constant, the equation (3.10) may be 


written 
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Inspection of last term in Eq. (3.12) and in view of [14] it seems ap- 


propriate to express it as 
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The magnetic field may be assumed to have the following form: 


Bed) =a ences (3.14) 
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Substituting (3.13) and (3.14) into (3.12) we obtain 
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On expansion of (3.15) the following set of differential equations are 


obtained 
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From the boundary conditions (3.9) and the equation (3.11) the following 


set of conditions for function f are obtained. 


f(0) = f5(0) = F5(1) = 0, f(1) = ¢ (3.21) 
f,(0) = (1) = F111) = 0, FH(0) = (3.22) 

and 
f(0) = f,(1) = f4(0) = 4(1) = 0 (3.23) 
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External Circuit Conditions: 

The unknown Ee in equation (3.13) may be obtained by considering 
the external circuit conditions of the bearing. If the bearing surfaces 
are insulators the total current flowing out of the bearing may be 


expressed as 


be lg 
I = | J, dxdy = o f ES dxdy + o | f U BY dydx (3.24) 
A A 0 O 


In terms of nondimensional parameters it gives 


C 
Tee = E, | dxdy + 0 | (2 + C,40,6' +...) dx (3.25) 
LUYuc h 


For the open circuit case I = 0 and equation (3.25) reduces to 


] 
C 
-E, | dxdy = 0 | (2+ cyteye! +...) dk (3.26) 
A ) 


Load Capacity: 
The bearing load capacity is found by integrating the pressure 
distribution. In terms of dimensionless variables the load capacity 


per unit width of the bearing is 


W= Le | mus (3.27) 


’ al 7 ’ , 
ota ah a 5 


a 
ex ve ite pia uy in | on pa 
eerie are eer... 7 nk 
ies cea: Sana a a 


aca at i807 a re es ies 


od 
tle): So 


sey abv | |e Taal haere |=" 
is 6 68 So 
(€5-&) wb... 04.0.4 E] er om =" 
a. }- baal ere 
| ae 
08 epauber (2S.£) rte be 0+ T ne ss 
- Bae - cy 
‘@S8) ° xb ae 0 ee an [erway i 3 e . (LU | 7 


aiuzeoq st pnizenpstat yd bauct at Gace, anid ait. _ 
wWiosyes baof siz arene" cata, Sele emvat fl : | 


(N88) | % | je 


: AU to i 


21 
3.4 Method of Solution 

The equations (3.16) through (3.20) are ordinary differential 
equations of third order and linear except equation (3.16). These 
equations and the boundary condi ate (3.21) through (3.23) are similar 
to the set of equations solved by Rodkiewicz and Anwar [14], except 
there are additional terms due to the magnetic effects. However, 
the same numerical methods can be successfully adopted. 

In order to obtain the solution to the problem, the dis- 
tribution of the applied magnetic field has to be assumed. In the 
present analysis, we shall consider three cases, Fig. 3.1, namely: 

(A) uniform field, 

(B) linearly increasing field, and 

(C) field proportional to (1 - Et ; 
Knowing the distribution of the magnetic field we can find the C's from 
equation (3.14), (see Appendix A) which are then substituted in 
equations (3.16) through (3.20) and (3.26). The solution of these 
equations give the constants D's. The pressure distribution is then 


obtained from expression (3.13). The load capacity of the bearing is 


given by equation (3.27). 


3.5 Discussion and Conclusions 

The pressure distribution curves for uniform and linearly 
increasing field are shown in Fig. 3.2. The M=0 curve, the case of no 
magnetic field, is also shown for comparison. It is noted that the 


position of maximum pressure changes slightly due to the presence of 
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magnetic field. It is evident from the figure that there is an in- 
crease in pressure distribution when the constant magnetic field is 
applied. Further increase is obtained when the applied field is made 
nonuniform. Consequently, the integration of the various pressure 
distribution curves results in the functional dependence of the load 
Capacity on the magnetic field distribution and the Hartmann number. 
These results are presented in Figs. 3.3 and 3.4. 

Figure 3.3 presents the relationship between the load 
capacity and the Hartmann number for neglected and included inertia 
terms. It is noted that the increase in load capacity due to the 
inertia effects decreases with the increase of Hartmann number. The 
no inertia terms results compare favourably with the low Hartmann 
number solution of Prakash [9], which extends only up to M= 1.2. The 
solution of the present work covers the range 0 <M < 5. 

In Fig. 3.4 the functional dependence of the ratio of MHD 
load capacity to the ordinary hydrodynamic load capacity is given. It 
is evident that the load capacity ratio for the nonuniform field is 
higher than for the uniform field for a given inlet Hartmann number. 
This increase depends upon the distribution function and the outlet 
Hartmann number. 

To show that the nonuniform field is more advantageous than 
the uniform field, we consider the case of uniformly increasing field 
of Hartmann number M. For such an applied field the average Hartmann 
number ratio MZ/M, which will be different than one, may be estimated. 


It is found that the load capacity ratio for this averaged Hartmann 
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ra) 
number, applied throughout the bearing, is smaller than the load capa- 
city ratio for the original linearly increasing field. 
For case C the distribution function is given by C(1 - 7 


In the present analysis, Cy was made equal to 1 giving M 2M. 


outlet — 
The corresponding load capacity curve is also shown in Fig. 3.4. 

From the load capacity curves it is noted that the curves 
for case B were terminated at lower Hartmann numbers than for the 
other cases. This was due to the convergence difficulties. 

The confergence of the power series (3.11), without the 
magnetic effects, was discussed by Rodkiewicz and Anwar [14]. It is 
noted that when the magnetic effects are included the convergence de- 
pends upon Hartmann number and distribution of the magnetic field. 

For case A and C, it was found that for Hartmann number up 
to 5, the series can be approximated by the first eight terms of 
(3.11). The contribution of the 9th and the higher terms were found 
negligible. 

It was noted that the convergence for case B, depended upon 
the gradient of the applied field and the Hartmann number. With the 
increase of these quantities it was found that more terms had to be 
taken into consideration in order to obtain the required convergence. 
For example, for Mp = 2 and M = 1.5, it was necessary to include up 
to fifteen terms to obtain reasonably accurate solution. 

Since the inertia effects on load capacity are noted to 
be small it is proposed that the inertia terms may be neglected. In 


Chapter IV the analysis will be considered without the inertia terms. 
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FIG. 3.3 INERTIA EFFECTS ON LOAD CAPACITY 
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CHAPTER IV 
EFFECTS OF MAGNETIC FIELDS IN SLIDER 
BEARINGS OF ARBITRARY PROFILE* 


4.1 Introduction 

In the previous chapter, the solution of an MHD bearing with 
nonuniformly applied magnetic fields was presented. In the analysis 
the inertia terms were retained and the resulting nonlinear partial 
differential equation was solved by a power series method. The solu- 
tion was limited to low Hartmann numbers due to the slow convergence 
of the proposed power series for the stream function. In the present 
chapter the inertia terms are neglected and the solution for the re- 


Sulting equations are presented. 


4.2 The Governing Equations and the Associated Boundary Conditions 


Since the inertia terms are neglected the governing equations 


(3.1) through (3.3) yield 


2 
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*n portion of this Chapter has been accepted for publication and will 
appear in the Journal of Lubrication Technology, ASME. July, 1972. 
And another portion has been accepted. for publication in Wear, an 
International Journal on the Science and Technology of Friction, 
Lubrication and Wear. (Forthcoming) 
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(f..8) “a (8 w+ a)o- os + ae = 0 | 7 


eae Belts 42 EY. (4.3) 


with the following associated boundary conditions: 


u(x,0) = U 


4.3 Equations in Dimensionless Form 


The above equations may be nondimensionalized by letting: 


B 
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The nondimensionalization will transfer the original fluid field into 


the square field as shown in Fig. 4.1. Equation (4.1) and (4.2) 


become: 


2— aes 
ou Te yuow o oe Tein oe mee ithe) 
—=z - h M B u-h MB =n; —- = 0 (4.6) 
aye y Jag? dx 
Muy LV. 9 (4.7) 
ox h ay 


The boundary conditions: 


u(x,0) = 1 
u(x,1) = v(x,0) = v(x,1) = 0 (4.8) 
AO) = culon = 
Defining the stream function such that 
Tees poe pee (4.9) 
oy OX 
Equation (4.6) becomes 
3— ys a open t = 
ov pew Be eee EL - W = 0 (4.10) 
oy y oy y dx 


4 Te oe Ps, ava : ip 7 . | | 
Pe A aa ; 
- am ma . aidest in 

| mY i. pene i py | 


a sis 
an 


(a8) eBig paar 


t= (0.%)0 


(8.8) 0.= (TRV = Ox = (120 


(@.0) 


With the following boundary conditions on jp: 
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External circuit conditions: 

To determine the unknown Ey in equation (4.10) it will be 
necessary to consider the external circuit of the bearing system. It 
is noted that when one of the bearing surfaces is a conductor it is a 


short circuited case which gives 
E_ = 0 (4.12) 


when the bearing surfaces are insulated and the end plates are perfect 


conductors we have from equation (4.3) 
Is | o(E, + u By) dx dy (4.13) 


or interms of dimensionless parameters 


iene 1 
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For the open circuit case I = 0 and equation (4.14) reduces to 


1 


=< [ [ (x) dx dy = WO] BY ok (4.15) 
Oo .O ce) 


Load Capacity: 
The load capacity may be found by integrating the pressure 
over the length of the slider. In terms of dimensionless variables, 


the load capacity per unit width of the bearing is 
~ Oae —) os 
W=—->= | Deax (4.16) 


Frictional force: 


The frictional force at the slider is given by 
Fee | ah: (4.17) 


In terms of dimensionless variables the frictional force per unit 


width of bearing is 


1 

Fh — 

(Prcmggeallls poate Akal be » ng (4.18) 
yu U Te 
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The friction factor Ce is given by 


(4.19) 
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4.4 Method of Solution 
The equation (4.10) is a linear differential equation of third 
order in which the coefficients are function of x. The equation is 
difficult to integrate directly in its present form. However, the 
solution to the problem can be obtained by considering the differential 
equation (4.10) at each section along the length of the bearing. At 


any particular section the general solution may be written as 


= vo y -va y oo 
We o, + fore + Cie Pa (4.20) 
where 
o(x)= fe Me BY (4.21) 
and 
a(x)= AX MB. EL + ne 22 (4.22) 
Aged dx 


using the boundary conditions (4.11) in (4.20) and eliminating the 


integration constants we obtain 


§ o3/¢ 1-02" - ha(lte 


2Vo (4.23) 


B(x)= - 


h(2+2e ie ence 


The foregoing values of 6(x) may be calculated if we know h(x) and Q. 
The values of h(x) are prescribed by the bearing profile while the 
value of Q is assumed in order to start the computations. Rearranging 


equation (4.22) we have 
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dp _ B(x) 7 OE 
r= Sts -MB E (4.24) 
dx ig MA Vc 


This equation (4.24) may be integrated numerically as the right-hand 
side is known at each point of the grid shown in Fig. 4.1. The nu- 
merically obtained value of pressure at the outlet is compared with 
the boundary condition (4.8) and if these two values agree within the 
specified limit then we have obtained the solution to the problem. In 
case of disagreement, a new assumption of Q is made and the above 
procedure is repeated. 

Once the pressure distribution in the bearing is obtained, 
the bearing load capacity and frictional force are obtained by nu- 


merical integration of equations (4.16) and (4.18) repsectively. 


4.5 Results and Conclusions 
The foregoing analysis and the method of solution may be 

applied to obtain the load capacity of a slider bearing of an arbitrary 
geometry where the distribution of transversely applied magnetic field 
may take any form of distribution along the length of bearing. How- 
ever, in this presentation a plane tappered type slider bearing is 
considered in the presence of an applied magnetic field which may 
have the following distribution (See Fig. 4.2): 

(A) uniform 

(B) linearly increasing 

(C) step-type. 


The results of these three cases are summarized below. 
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(A) Uniformly Applied Magnetic Field 

As pointed out in [3,4], for uniformly applied magnetic fields, 
the greatest increase in load capacity occurs if the bearing surfaces 
are perfect electrical insulators. Thus, the same properties are as- 
sumed here. For any Hartmann number the load capacity may be obtained 
by fixing the ratio of film thickness at inlet to outlet hy. However , 
it is noted that for a given Hartmann number there is an optimum value 
of is for which the load capacity obtained is the greatest. Hence, 
the calculations were carried out using various values of hye for any 
Hartmann number, until the optimum hy. was found. The values of the 
load capacity and the optimum h,. are plotted against M in Fig. 4.3. The 
corresponding plot of the flow rate through the bearing Q vs M is 
shown in Fig. 4.4. It can be seen that the optimum h. increases with 
the increase of M at low Hartmann numbers and reaches its maximum 
at M= 3. Any further increase in M results in the decrease of optimum 
hh. The optimized load capacity increases with the increase of M. 

In order to show the influence of optimum ne and conductivity 
of the bearing surfaces the results of Shukla [8] are also shown in 
Figs. 4.3 and 4.4. It is noted that the load capacity obtained in the 
present analysis is slightly less than obtained by Shukla at low M, 
but are higher for large Hartmann numbers. Therefore Shukla's arrange- 
ment may be considered slightly advantageous. At higher Hartmann 
numbers, however, the present configuration with insulated surfaces 
gives better results when the film thickness is optimized. Thus the 
stator surface need not possess variable conductivity (which might 


be difficult to achieve from the practical point of view) in order to 
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obtain the greatest increase in load capacity. Figure 4.4 shows that 
for M > 2,Q decreases with the increase of M in both analyses. However, 
there is a considerable difference between the two values at any Hart- 
mann number although the minimum film thickness is the same. This 
difference may be explained by observing that both bearings form two 
different systems for the induced currents. As a result the effects 

of body force on flow are different for each of the cases. For Shukla's 
arrangement the effects of this force is to retard the flow near the 
outlet of the bearing. In the present system an accelerating force 

is generated near the inlet and retarding force near the outlet. This 
arrangement reduces Q, but the overall effects are not as dominant 

as in Shukla's analysis. 

In Fig. 4.5 the relation of the frictional force and the 
friction factor, to the Hartmann numbers is indicated for the same 
optimum conditions as Fig. 4.4. It can be seen that the frictional 
force increases with the increase of M while the friction factor de- 
creases. 

(B) Linearly Increasing Applied Magnetic Field 

Increasing applied magnetic field, along the length of the 
bearing, may take any form i.e., linearly increasing, parabolic dis- 
tribution etc. In the present analysis the calculations are presented 
for the linearly increasing field in order to show the application of 
the present method to nonuniform magnetic fields. The slider bearing 
considered is assumed to have insulated surfaces. The distribution of 
the linearly applied magnetic field is assumed such that the intensity 


of the magnetic field at outlet is twice that at the inlet of the 
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bearing. This gives the ratio of M M, = 2. In obtain- 


outlet! “iniet ray | 

ing the solution for any Hartmann number the ratio i is optimized, 

as described in Case A, in order to obtain the greatest load capacity. 
A plot of the optimized hy. and the load capacity versus M is presented 
in Fig. 4.6. It is observed that h, decreases with the increase of 
M(M>1) while the load capacity increases. The load capacity obtained 
for the uniformly applied field, Case A, is also shown in the same 
figure. From Chapter III and the present results, it is noted that 
nonuniform applied magnetic fields give higher load capacity than 
comparable uniform magnetic fields. 

A plot of frictional force and friction factor against M 
is shown in Fig. 4.7. It can be seen that the trends of these curves 
are similar to these noted for Case A. 

The difficulty in obtaining the solution for high Hartmann 
numbers which was reported in Chapter III has been eliminated with 
the present method by ignoring the inertia terms. 

(C) Step-type Applied Magnetic Field 

For a step-type case the magnetic field of uniform strength 
is applied in the downstream part of the bearing. The stator of the 
bearing is assumed of conducting material while the bearing is an in- 
sulator. It is noted [7] that a bearing of such an arrangement 
generates load carrying capacity even when the two surfaces are parallel. 
The load capacity depends upon the Hartmann number, applied magnetic 


field step location S, and the ratio of the film thickness at inlet 


to outlet, hy However, it is noted that for a given Hartmann number 
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there is an optimum value of ie and S for which the load capacity is the 
highest. Hence the calculations were carried out using various values 
of h, and S, for any Hartmann number, until the optimum h, and S were 
found. The optimized ne and the corresponding load capacity are 
plotted against M in Fig. 4.8. It can be seen that the optimum hy. 
decreases with the increase of M. It becomes equal to one when M = 4. 
This is the case of the parallel slider bearing. When M > 4 the bear- 
ing film becomes divergent. The optimum step location is fairly con- 
stant (S = .65) when M < 3. Further increase in M decreases the step 
location until M = 5, and thereafter it becomes practically constant 
again (S = 0.55). The results of Kuzma [7] are also shown in the 
same figure. It can be seen that when the ratio hy. is optimized the 
maximum load capacity is higher than for the parallel slider bearing. 
This difference becomes quite significant at large M. It is interesting 
to note that when M = 4 and S = 0.6, the parallel slider bearing gives 
the highest load capacity. 

Figure 4.9 shows a plot of the frictional force and the 
friction factor, for maximum load conditions, as function of M. It 
can be seen that the frictional force increases with the increase of 
Hartmann number while the friction factor decreases. For comparison 
the plot indicates the results of Kuzma [7]. It is noted that the 
values of frictional force obtained, for M > 5, are higher than those 
of a parallel slider bearing. However, the friction factor values 


are reasonably in agreement. 
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Conclusions: 

A simple method of calculating the load capacity for an MHD 
Slider bearing has been described in this chapter. The transversely 
applied magnetic field may have any form of distribution along the 
length of the bearing. The method may be conveniently applied to 
slider bearings of a general profile where the bearing surfaces may 
be insulators, conductors, or possess variable conductivity in x- 
direction. Foom the results the following conclusions are drawn: 

(1) For uniformly applied magnetic fields the load capacity 
obtained is the highest if the ratio of film thickness at the inlet 
to that at the outlet is optimized and the bearing surfaces are insu- 
lated. 

(2) For step-type applied magnetic fields one obtains an increase 
in load capacity if the stator is a conductor. However, to obtain 
the maximum increase the ratio of film thickness at the inlet to 
that at the outlet has to be optimized. 

(3) For linearly increasing applied magnetic fields the highest 
load capacity is obtained if the ratio of film thickness at the inlet 
to that at the outlet is optimized and the bearing surfaces are insu- 
lated. The results indicate that the nonuniform magnetic fields give 


higher load capacity than the comparable uniform magnetic fields. 
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CHAPTER V 
NONUNIFORM ELECTRIC FIELD EFFECTS 
IN MHD SLIDER BEARINGS 


5.1 Introduction 

The review of literature and the results presented in the 
previous chapters indicate that at smal] Hartmann numbers only a slight 
increase in load capacity is obtained under open circuit condition. 
However, a Significant increase may be obtained when Hartmann numbers 
are very large. Thus, to substantially increase the load capacity a 
very high strength magnetic field would be required. Consequently, this 
would necessitate a very large magnet which would be impractical. 

In order to improve the load capacity at low M the supply of 
electrical power to the bearing, by connecting the side electrodes to 
the external source, has been proposed [5,6,9,10]. This arrangement 
has substantially improved the bearing load capacity even at very low 
Hartmann numbers. 

However, in these analyses it was assumed that since the 
electrodes were perfect conductors, and the bearings were assumed 
infinitely wide, then from equations (2.3) E was essentially a constant. 
Therefore, any variation of E. in the bearing and neat the edges was 


neglected. 


It was concluded in Chapters III and IV, that the nonuniformly 
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applied magnetic fields improve the average load bearing capacity. This 
intuitively suggests that if the applied electric field from a certain 
part of the bearing is removed, then further improvement in the load 
Capacity can be expected. This may be achieved by segmenting the side 
plates into the conductor and insulated parts. The external power is 
then only supplied to the conductor part. 

The purpose of this Chapter is to investigate the effects 


on load capacity of such a bearing configuration. 


5.2 The Bearing Configuration 


The finite width parallel plate bearing with a uniform mag- 
netic field applied perpendicularly to the bearing surfaces is shown 
in Fig. 5.1. The bearing surfaces are made of insulating material. 
Each of the side plates is segmented into two parts; the electrode, and 
the insulated. The electrode part is a perfect conductor, while the 
latter part is an insulator. The electrode part may be connected to 
the external generator such that the current can be supplied to the 


bearing. 


5.3 The Governing Equations and the Associated Boundary Conditions 


The general equations developed and simplified in Chapter I1, 
remain valid for the present model. However, when the external potential 
is applied to the bearing system shown in Fig. 5.1, it becomes necessary 
to consider the Maxwell's curl equations as well. Thus, to obtain the 
solution all these equations have to be solved simultaneously subject 


to the appropriate boundary conditions on the fluid and electric fields. 
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This indeed, is a difficult task and at this stage we make some further 
assumptions to put the problem in a manageable form. These assumptions 
are: 
(1) It may be recalled that the coupling of the fluid equations 
and Maxwell's equations are through Ohm's Law, stated in Chapter II 


as: 


J = o(E + V x B) (2.4) 


In the present analysis consideration is essentially on the MHD where 
the applied magnetic field is not high, and the velocity of the slider 
is not large. Therefore, it would be reasonable to assume that the 
induced voltage, V x B is much smaller than the applied voltage. Hence, 
for such a case the second term may be neglected compared to the first 
term in equation (2.4), and the equation may be written in the following 


form: 


J = o(E) (5.1) 


By making the above assumptions, the fluid equations have been uncoupled 
from Maxwell's equations. Therefore, each set of equations can be 
solved independently with the appropriate boundary conditions. 

(2) For the bearing under consideration E and J may be assumed 
to be functions of x and Z. 


(3) As indicated in the results of Chapter III, the inertia terms 


are neglected. 
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From assumptions (1) and (2) and equation (2.3) the follow- 


ing results: 


V.E=0 (no) 


Defining Eee Vo, and combining this definition with the above equation, 


the Laplace's equation for the two dimensional domain (x,z) is obtained. 


Ey a0. 4G (5.3) 


The solution of this equation with the appropriate boundary 
conditions will give the ¢ distribution, and hence the electric field 
in the bearing. 

In order to solve equation (5.3) it is necessary to specify 
the boundary conditions on the four sides of the domain xz. From 
Fig. 5.1, it is apparent that the electric potential is applied to 
the electrode parts of the side plates, while the remaining parts 
of the side plates are insulators. The boundary conditions, however, 
at inlet and outlet of the bearing are difficult to prescribe due to the 
continuous flow of the lubricant. Since the slider bearing essentially 
represents a portion of the thrust bearing arrangement, it may be noted 
that at the end of each thrust pad the groves of small width for entry 
and discharge of the lubricant are provided. It may be pointed out 
that this arrangement can be modified so that near the inlet and the 


outlet the insulated walls may be provided. However, in the present 
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analysis, the arrangement is approximated by assuming that the insulated 
walls are located at the inlet and outlet of the bearing. The boundary 


conditions for such an arrangement may be written as 


82 (x, +2) = 0 when x>L, (5.4) 
Dine 
o(x, #5) = +o, when x <L, 


Under the assumptions (1) through (3), the governing equations 


(2.6) through (2.11) now take the following form 
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From equations (5.5), (5.6) and (5.7) when combined with 


the boundary conditions (5.10), the magnetohydrodynamic form of 


a} 


Reynold's equation may be obtained by employing the following procedure. 


Since Jy and J, are assumed to be functions of x and z only, 


equations (5.5) and (5.6) may be integrated directly to yield 
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where C's are constantsof integration. When these are evaluated with 


the use of boundary conditions (5.10) then 
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The x-component of the fluid flow rate per unit width, Qy 


can now be expressed as 


h 


QU f udy (5.16) 
(@) 


Combining (5.14) and (5.16), the following is obtained 


3 
__ hap Uh 
Whe? Span toys oes eg) ar (Sal) 


Similarly, the z-component of fluid flow rate per unit 


length is expressed as 
Q_ = | wdy (5.18) 
which, when combined with equation (5.15) yields 
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At this stage the continuity equation (5.7), may be inte- 
grated with respect to y. Using boundary conditions (5.10), the follow- 
ing is obtained. 

h h 


9 5 : 
x | u dy + val w dy = 0 (S220) 
0 0 


From equations (5.16) and (5.18) the above equation becomes 
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Substituting equations (5.17) and (5.19) into equation (5.21) yields 
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When h is not a function of x as is the present case, the above equation 


reduces to 
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It is interesting to note that the above pressure equation 
is independent of the viscosity of the lubricant. At this point it 
is convenient to put the governing equations into dimensionless form 
by using the characteristics quantities usually employed in MHD lubri- 


cation. These are: 
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In dimensionless form the potential equation (5.3) becomes 


8 (x, +5) = 0 when x > S 
dZ 
g(x, +5) 49, when x < S 


While the pressure equation (5.23) reduces to 
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Load Capacity: 


The integration of the pressure over the bearing area yields 


the load capacity. In terms of dimensionless variables it may be 


written as 


Frictional Force: 


The frictional force on the slider is given by 


ype 
eee f | u(t) dx dz 
ay y=0 
-D/2 0 


From equation (5.15) the above equation becomes 


D/2aeel 
r= | | Ly (B+ Bg) + HU] dx dz 
-D/2 0 


which when simplified gives 
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5.4 Solution to the Governing Equations 
5.4.1 General Approach 


The problem formulated in section 5.2, requires the solution 
of equations (5.25) and (5.27) with the corresponding boundary condi- 
tions (5.26) and (5.28) respectively. These equations are difficult 
to solve analytically and hence a numerical approach is adopted. 

In the present study, the method of finite-difference will 
be used [12,15-17]. The resulting equations will be solved using the 


point successive-overrelaxation method [17]. 


5.4.2 Finite-Difference Approximation 


Considering a 5-point molecule shown in Fig. 5.2, the partial 
differential equation (5.25) is transformed into a finite-difference 
form (see Appendix B). After rearranging, the equation for any point 


(i,j) in the grid may be written for Ax = Ay, in the following form 
(5.33) 
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Similarly using the 5-point molecule shown in Fig. 5.3, 
equation (5.27) for Ax = Ay = Ah, becomes (see Appendix C) 
D. b 


Meal 7 CMD 
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5.4.3 Method of Solution 


The finite-difference approximation at each of the grid points 


has essentially reduced the problem to solving a set of ordinary linear 


equations. 
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From Fig. 5.2, it is clear that due to symmetry about the 
middle line only one half of the domain may be considered for computations. 
(a) Solution for Electric Potential Equation 

Since the equation is uncoupled, the solution to this equation 
(5.25) with boundary conditions (5.26) is obtained first. 

Due to unusual boundary conditions on ¢ the half domain was 
further divided into two regions; close to the electrode and far away 
from it. In region I, close to the electrode it was divided into a 
finer mesh as compared to the far away region. This was primarily to 
overcome the difficulties in computations due to the steep gradient 
of » near the electrode plate and especially the corner. 

The application of equation (5.25) at each of the grid points 
gives a set of linear equations. These equations are then solved by 
the successive-overrelaxation procedure [17] by a digital computer, 

IBM 360/67. 
If the number of iterations performed is denoted by K, the 


procedure can be defined with: 


(5.35) 


where w iS a parameter known as a relaxation factor, the choice of which 
determines the rate of convergence. The iteration procedure is terminated 


when the following convergence condition is satisfied. 
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Once the potential distributing in the bearing is obtained, 


J. 


J be 


ee 3J,/dz and aJ,/ax may be approximated at each of the 
grid points by the finite-difference formulas given in Appendix D. The 
total current through the bearing for the applied potential %, is then 
obtained by carrying out the numerical integration of J, distribution 
at the center line. The area to be considered for the total current 
will be in the xy plane. 
(b) Solution for Pressure Equation 

Knowing the electric field and the current distribution in 
the bearing, the pressure equation (5.27) may be solved. A uniform 
grid size shown in Fig. 5.3, is used in calculations. The procedure 


of solving the finite-difference equations is the same as described 


above. Equations (5.35) and (5.36) may be written in terms of pressure 


as 
_K _(K-1) _(K-1) — _(K) EA K= 10 aK) 
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Once the solution to the pressure equation is obtained, the 
load capacity and the frictional force on the slider are found by inte- 
grating the equations (5.29) and (5.32) respectively. The method 
employed for integration was Simpson's Integration Rule. 

The optimum values of relaxation factor w, and the mesh 
size employed in the computations is given in Appendix E. The error 
involved in the numerical solution is also briefly discussed in this 


Appendix. 


5.5 Results and Conclusions 

As shown in Fig. 5.1, the electrical power is supplied to the 
bearing by connecting the side electrodes to the external source. For 
an applied potential to the electrodes, 5: the electrical current 
flowing through the bearing system will depend on the electrode length 
S. In order to study the effects of nonuniform electric field on the 
bearing load capacity, the current flowing through the bearing is kept 
constant, while the electrode length S is varied. The constant flow 
of current can be achieved by adjusting the value of the applied 
potential o,- The results of such a bearing arrangement with L/D = 1.0 
are discussed below for M = 0.1 to 1.0. 

The solution of equation (5.25) gives the potential distribu- 
tion in the bearing. The current density in the bearing is then ob- 
tained by using the formulas given in Appendix D. Knowing the current 
density distribution, the current streamlines are obtained by inte- 
grating the current density distribution. The potential lines for 


5 = 2.02 and S = 0.4 are shown in Fig. 5.4. For the same conditions, 
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Fig. 5.5 shows the current streamlines in the bearing. The tendency of 
concentration of electric currents near the corner of the electrode is 
exhibited in this figure. It was noted that with the increase of the 
electrode length, this concentration of current decreased. 

The solution of the pressure equation (5.27) presented in 
Fig. 5.6 indicates that there is a pressure buildup in the bearing 
when an electric potential is applied to the electrode parts of the 
side plates. The results further reveal that for a given hydrodynamic 
condition and for a fixed current flow through the bearing, the pressure 
distribution depends upon the electrode length S, and the Hartmann 
number. The pressure buildup in the parallel plate bearing with seg- 
mented side electrodes may be explained as follows. 

From the current streamlines diagram (Fig. 5.5) it can be 
seen that when the electric potential is applied to the side electrodes, 
there is a net current flow in the negative direction of z-axis. The 
interaction between the current passing through the fluid and the 
applied magnetic field results in an electromagnetic body force in 
the fluid film. In the present analysis, the important components of 
the electromagnetic force as shown in equation (5.5) and (5.6) are: 
J, x B. in x-direction, and Jy X BY in z-direction. The distribution 
of the body force in the bearing for uniformly applied magnetic field 
therefore will depend on the distribution of the current density. 
In order to illustrate the role of the electromagnetic body force on 
the flow more clearly, let us consider the situation at the center of 


the bearing. As the bearing is symmetrical about the middle line, Jy 
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will be zero at the center. Therefore, the only component of the body 
force at the center will be the x-component. If one visualizes the 
direction of this body force, it can be seen that the effects of the 
electromagnetic force is to accelerate the fluid between the bearing 
Surfaces in the same direction as the motion of the slider, thus re- 
sulting in an increase in pressure. However, if the electromagnetic 
force is distributed uniformly, then there is no pressure buildup in 
the parallel plate bearing. Mathematically, this can be easily seen 
from equation (5.27) and the boundary conditions (5.28). For a uniform 
electric field, the right hand side of the equation is zero, so that 
p = 0 throughout the bearing. However, for a nonuniform electric 
field, the right hand side is nonzero resulting in a nonzero pressure 
distribution. 

The physical reasoning of the pressure rise in the bearing 
may be explained if the MHD parallel plate bearing is visualized as 
consisting of two parts. For simplicity, let us assume that the re- 
sulting electric field is uniform and is only confined to the electrode 
part of the bearing. Then, the bearing may be divided into two distinct 
parts; the upper part where the fluid is being pumped by the electro- 
magnetic body force, and the lower part which essentially represents 
an ordinarily couette flow superimposed by a forced flow. In order 
that the couette part sustains the extra fluid which is being pumped 
through the bearing, a definite pressure drop will be required. There- 
fore, a positive pressure gradient develops in the pumping part of the 


bearing. However, this increase in pressure will depend on the con- 
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ditions in the lower part. The situation may be viewed as somewhat 
similar to a step-type hydrodynamic bearing. As will be seen later, 
that the pressure developed in the bearing becomes smaller, if the 
couette part is made smaller by increasing the length of the side 
electrodes. Finally, the pressure reduces to zero when the side 
electrodes are extended to the full length of the bearing. Then, the 
bearing configuration simply represents an MHD couette arrangement. 

Figure 5.6 shows the pressure distribution along the length 
of the bearing for I = 2.5, S = 0.1, and M= 1.0. The curves are 
plotted for the center of the bearing (z = 0), and z = 0.25 and 0.375. 
It is noted that the maximum pressure distribution occurs at the center 
of the bearing indicating the location of the maximum electromagnetic 
body force in x-direction in the bearing. The curves further show 
that the position of maximum pressure moves towards the inlet of the 
bearing as z increases from the centerline. The reason for this shift 
may be attributed to the concentration of current near the electrode. 

The load capacity of the bearing against M for various values 
of S is shown in Fig. 5.7. It is noted that for a fixed S, the load 
capacity increases with the increase of M. The effect of electrode 
length S on load capacity is clearly indicated in this figure. It 
can be seen that as S is increased, the pressure developed in the 
bearing becomes smaller. Finally, for S = 1.0, the pressure reduces 
to zero, and the bearing has no load carrying capacity. 

Figure 5.8 shows a plot of frictional force on the slider 


against M. It is noted that the frictional force increases with the 
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increase of M. However, it may be pointed out that the frictional force 
is independent of the electrode length S. Finally, the frictional 
force as shown in equation (5.32), depends on the total current through 
the bearing for a given Hartmann number. 

The foregoing results of MHD parallel plate slider bearing 
suggest that if the electrical power is supplied to the bearing to 
improve the load carrying capacity, it would be advantageous to apply 


this power near the inlet portion of the bearing. 
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CHAPTER VI 
CONCLUSIONS AND SUGGESTIONS FOR FUTURE INVESTIGATION 


6.1 General Results and Conclusions 

A theoretical study was made of the effects of electromagnetic 
interactions on the fluid film of a liquid-metal-lubricated slider 
bearing where the magnetic field was applied perpendicularly to the 
bearing surfaces. The analyses were carried out for an open circuit 
condition, and for the case where the electrical power from an external 
source was Supplied to the bearing. The main results are summarized 
below. 

For an open circuit condition the applied magnetic fields 
were assumed to have uniform, increasing, and step-type distribution. 
The effects of inertia on load capacity were analysed. The results 
indicate that the contribution of inertia terms decreases with the 
increase of Hartmann numbers, and it becomes negligible at high M. 
Therefore, the inertia terms were neglected and the solution was pre- 
sented for uniform, linearly increasing and step-type applied magnetic 
fields. It is noted that for each type of magnetic fields, the maxi- 
mum load capacity of bearings is obtained when the ratio of film 
thickness at the inlet to that at the outlet is optimized. The bear- 
ing surfaces have to be insulators for uniform and increasing magnetic 


fields, while the stator have to be a conductor for a step-type magnetic 
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field. The results also indicate that nonuniform magnetic fields i.e., 
increasing type, give higher load capacity than the comparable uniform 
magnetic field. However, the improvement obtained in load capacity by 
applying a nonuniform magnetic field is only moderate. For step-type 
applied magnetic field the load capacity obtained is less as compared 
to the uniformly applied magnetic field. It is interesting to note 
that for such a case, when M > 4 the maximum load capacity is obtained 
when the fluid film is divergent. The method of solution developed 
in Chapter IV may be employed in analyzing bearings having complicated 
profiles, e.g., multistepped and composite bearings. 

For the case of externally supplied electric current, a 
parallel plate slider bearing was analyzed. The magnetic field was 
applied perpendicularly to the bearing surfaces. The segmented side 
electrodes were connected to the electrical power source, thus giving 
a nonuniform electric field distribution in the bearing. The results 
indicate that a parallel plate bearing with such an arrangement have 
load carrying capacity. The results suggest that if the electrical 
power is supplied to the bearing to improve the load capacity, it 
would be advantageous to apply the power near the inlet portion of 


the bearing. 


6.2 Suggestions for Future Investigation 


In the present work the effects of nonuniformly applied magnetic 
and electric fields were analyzed for a simple slider bearing. The 


analyses may be extended for bearings with complicated profiles. 
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For nonuniformly applied electric fields considered in 
Chapter V, the fluid equations were uncoupled from Maxwell's equations 
by making certain assumptions. It is suggested that these equations 
may be solved in the coupled form and the results compared with those 
presented here. 

Another possible aspect of further investigation in magneto- 
hydrodynamic lubrication is attempting to design a bearing having 
automatic electrical adjustments in order to compensate for transient 
loads. 

The present investigation and the results of other authors 
as discussed in review of literature, have shown the possibility of MHD 
slider bearing. In order to develop it into a practical device, it is 


suggested that a series of experiments should be carried out. 
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APPENDIX A 
MAGNETIC FIELD DISTRIBUTION 


The expression for the magnetic field, equation (3.14), may 


be written 
ig si + C56 + C8° + (Awl) 
y § ] 2 3 eee e 
For Case A 
C, = Lege Co = Cy = C, mes 0 (Ase) 
For Case B 
C, ca Seman C, = - 2S 
(A.3) 
Co = C4 = Cy eee) 
dB. 
where Sa = —t 
dx 
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APPENDIX B 
ELECTRIC POTENTIAL EQUATION IN 
FINITE-DIFFERENCE FORM 


The potential equation (5.25) may be written for convenience 


as 


CER 7a) ee (5.25) 
axe D” gaze 


Considering the five-point molecule shown in Fig. 5.2, the following 


may be written [19] 
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Substituting (B.1) and (B.2) into equation (5.25), the follow- 


ing is obtained 
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Az = Ah, the above equation becomes 
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APPENDIX C 
PRESSURE EQUATION IN FINITE-DIFFERENCE FORM 


Rewriting the pressure equation (5.27) from Chapter V as 


ne mere’ 
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Substituting (C.1) and (C.2) into equation (5.27) yields 
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Az = Ah, the above equation may be written as 
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APPENDIX D 
FINITE-DIFFERENCE FORMULAS TO OBTAIN ELECTRIC FIELD 
AND CURRENT DENSITY DISTRIBUTION IN THE BEARING 


Using the 5-point molecule shown in Fig. D.1 the following 


finite-difference formulas may be written: 
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FIG. D.L 5-POINT MOLECULE 
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Knowing a. and u, distribution in the bearing, (30/82) 
and (J, /0x) may be approximated at each of the grid point by using 
the same five-point molecule. The formulas may be easily written 
from the above expressions by replacing $'s with appropriate values 


of current densities. 
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APPENDIX E 
CONVERGENCE PARAMETER, MESH SIZE AND 
ERROR INVOLVED IN COMPUTATIONS 


E.1 Relaxation Factor For Point Successive-Overrelaxation Method 
For point successive overrelaxation method, the optimum 
values of relaxation factor w, are discussed in [17,18]. The important 


relation for linear elliptic type partial differential equations given 


in [17] are 
(i 
“optimum | craw cn ae (E.1) 
ll eee alles ee ul 
where A= x (cos = + cos b) : 


a/2 = Mh, b = Nh and h = mesh size. M and N are the number of divisions 
in x-direction and in z-direction respectively. 

In solving the potential equation (5.25) and the pressure 
equation (5.27), the optimum w for the grid size employed in computa- 


tions, was calculated from the above relation. 


E.2 Error and Mesh Size 
The error is inherent in the numerical solution. The errors 
due to finite difference approximations of equations (5.25) and (5.27) 


are of 0(A“h). These errors can be reduced by increasing the number 
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of divisions. However, the computing time increases considerably as 
the mesh size used is further reduced. In the present work, the follow- 
ing mesh size was employed which was in consistence with the tolerance 
specified in equations (5.36) and (5.38). 
For potential equation the mesh size for Region I and II 

was .0125 and 0.025 respectively, while for pressure equation it was 
nOZ5p 

Since the computations were carried out using the double 


precision numbers, the round off errors should be insignificant. 
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APPENDIX F 
NUMERICAL EXAMPLE FOR EXTERNALLY 
SUPPLIED ELECTRIC POWER 


The results of nonuniformly applied electric field are pre- 
sented for parallel plate slider bearing in Chapter V. To heip visua- 
1ize the load that can be obtained and the electrical power required to 
maintain this load, a numerical example will be considered. Typical 


numerical values of the parameter used are (RMKS units): 


h = 107¢m 

L = 0.1m 

D = 0.1m 

U = 10 m/sec 


Using mercury as the lubricant, the following properties of 


mercury are taken from reference [20]. 


u = 1.55 x 1073 Newton 5ec/m< 


1.07 x 10° mho/m 


Q 
iT} 


For S = 0.1, M= 0.5, and I = 2.5, the bearing load capacity 
Was obtained from Fig. 5.7, is 0.215. From these dimensionless numbers, 


the important physical quantities are calculated. These are: 
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APPENDIX G 
COMPUTER PROGRAMS 


‘ae ANT 
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2 ¥1GWS99A 
eMAHOORY AITUIMOT 


CeRKK 
& 
CKKK* 
Cc 

C ** 
Ce X 


OR gO tok ak ok a dak kb dak ak ok ook 


* IN TRE SIXTH COLUMN INDICATES 
CONTINUATICN FROM PREVIGUS LINE 


FE RK Be oh ae Oi a A a OK a Fe a he SK ie aie a OK a OK a a oe OK a a ak 


FORTRAN PROGRAM FOR CHAPTER 3 
THIS PROGRAM CALCULATES PRESSURE DISTRIBUTION IN 
MHD BEARING 
IN THE ANLYSIS THE INERTIA TERMS ARE RETAINED. THE 
SOLUTICN IS OBTAINED BY SOLVING A SERIES OF 3RD ORDER 
DIFFe EQUATIOGNS.~IN THIS PROGRAM ONLY SEVEN EQUATIONS 
ARE SOLVED 
SUBROUTINES FRCMW IBM-SYSTEM/360 SCIENTIFIC SUBROUTINE 
PACKAGE ARE USED IN THIS PROGRAM 
DGUBLE PRECISION F0(55)-2F01(55).,F02(55) 2F 03155) sF1(55) 
x 9F11(55)- 
1F121(55}3.F13(55)s REsHARTsCOsC1sC2sC3sALPHA(5),B015) 3G 
- AMA(5).CB> 
2PRMT (5) oe Y(3) saDERY(3) AUX (804) s A155) 38(55)2C(55).D0(55)— 
5 £(55}> 
3 2Z(55)5R155)4.CC(55)5AF1(5)5CF(5),.21(55).81(5).88s8215) 
% 2f2(55)>. 
4F21(55).2F22(55).F23(55).221(55) »Z11155).212(55),2Z213(55) 
* 2B3(5)6 
5F3(55)0F31(55),F32155).F 33155) »DEL( 21) +PGRAD{I21)5SsZPi{ 
x 21) sB4(5)- 
685(5) .F4(55) .F41(55).F42(55) sF43(55). F5(55)2F51155) sF 
5 52(55)>. 
7FS3(055),;CUR eEE(21) »386(5).87(5).88(5).F6(55)sF61(55)> 
*K FEZ20S5 }s 
8F63(55) 
EXTERNAL FCTs,OUTP 
COMMON FOeFO1sFO2eF03eRE 2HARTsCBsC02TIC 
COMMON /FIRST/SIL 
CALCULATION FOR 3RD ORDER DIFF e EQNS. 


SOLUTION FOR EGe NOs 1 1S GETAINED USING FOURTH CRDER 
RUNGE KUTTA METHOD 

CALCULATION FOR EQe NOol 

FILM THICKENESS RATIO = 220 


cts) 


Od HTX] ‘T UL * 
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efide € }s f * 
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CK 
& 
a2 
C 

27 
C#EKEX 
Cc 
Cxxx 

31 

45 

32 

33 


DG 23 MM=1.42 


REAC(5:25) ALPHA{MM) »sBO(MM), GAMA{MM) 


FORMAT {(3F15-9) 
RE=1.20 


THE APPLIED MAGNETIC FIELD IS SPECIFIED BY THE VALUES 


OF COsCisC2sC3rse0e 


WHEN CO=1.20 AND C1l=C2=C3=020 »sTHE APPLIED MAGNETIC 


FIELD IS PROPORTIONAL TO (1-X/2)**-1 


C0O=1.20 

C1i=0.0 

C2=0.0 

C3=020 

CB=BO(MM) 

BETA=ALPHA(MM) 

HART= GAWMA(MM) 
WRITE(6227) CBsBETAsHART 


FORMAT(70%,.5X%s *CB="sF1509% "BETA=* sF1l5e9s 


x 026) 
PRMT(1)=0.20D0 
PRMT(2)=1.20 
PRMT{(3)=2025 
PRMT(4)=020001 
¥{€1)=0.0D0 
¥(2)=0.20D0 
Y¥(3)= BETA 
DERY(1)=225D0 
DERY(2)=225D0 
DERY{13)=0.25D0 
NDIM=3 
iIc=0 


"HART="',F1i 


CALL DRKGS(PRMTs Y¥2DERY »NDIM, ITHLF »FCTsOUTPsAUX) 
EQUATIONS NOs 2 THROUGH 7 ARE PUT IN FINITE DIFF .FORM 
AND THE RESULTING RECURRANCE RELATION IS THEN SOLVED 


CALCULATION FOR EQe NOe 2 
Ti=41 

TliI=1t-1 

H=0 2025 

READ(5531) Bi( MM) 
FORMAT{(F 1026) 

BB=Bi1(MM) 

KK=0 

IL=1 

WRITE(6s32) 8B 
FORMAT(40'%,5Xs *Bl="s F1026) 
KK=KK+1 

WRITE(6333) 


FORMAT(70% .9Xs "A%s20XKs*B* sy 215X%0*C*% og 15X%9* D8 915X%o SE* ) 


DO 34 N=3,II1!I 
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47 


44 


43 


53 


48 


49 


50 


54 
46 
CO 


75 
62 


ACNI= ©1240 05*(CH**2) #( REXFO1 (NI— CHART **2 )*( CO*%2))) 


BIN)= (CH*¥*3)*F02(N) *RE-3.) 


CIN) = (30-0 25*(H**2) *(RE*FO1(NI-CHART**2) *(CO**2))) 


DIN) =—-1.-0 
E{N)= (H**3) *BB 
CONTINUE 


CALL SOLU(AsB2CsDsEsZZ) 
DO 47 N=2sII 
FI(N)=ZZ(N) 

CONTINUE 

F1(1)=0.0 
APFICKK)=FICIT) 

CF(KK)= BB 
TFCDABS(FITLTI))-« LTe 12E&-6) GO TO 43 
IF(KKeEGQe1) GO TQ 44 
TFC KK.EQ@e3) GO TQ 23 
DXX= AF1(2)-AF1(1) 
DYY=CF(2)-CF(1) 
BB=CF(1)—-AF1(01)*DYY/DXX 
GO TO 45 

BS=88+2055 

GO TO 45 

WRITE(6,53 ) 


FORMAT ©€2 0% .9X%s4F1%315X0e"F11% 21SXs'Fi2*» 


Bit(MM)=BB 

NDM=41i1 

CALL DDETS (HsesFil »sZsNDMSIER) 
DO 48 N=1,I1I1 

F1I1(N}=Z(N) 

CALL DDETS (HsF1ll»Zs2NDM,IER) 
DG 49 N=1,I1if 

Fi2(N)= ZINN) 

CALL DDETS5 (H»sF12+ZsNDMs IER) 
DO 50 N=1,11 

F13(N)=Z{N) 

DO 46 N=i1,I11I1 
WRITE(6,54)F1(N) sFI1 (Nd 2F12IN) sFI3(N) 
FORMAT( 40% »5X%s4F 1526) 

CONT INUVE 

CALCULATION FOR EQ. NGe 3 
IL=IL+4i1 

READ(5231) B2(MM) 

BB=B2(MM) 

KK=0 

WRITE(626€2) 8B 
FORMAT(* 0% »5Xs "B2="»s F106) 
KK=KK+1 

WRITE(6433) 
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DO 64 N=3,II1I 
ACN 9=€ 1 0— 0 S¥*¥(H*EK2 ) KL HART *#*¥2 9) *(CO**2) ) 
B(NJ= (H**3)*2*RE*XFO2(N)I-320 
CIN) =(324005%(HE2) *(HART** 2) *¥({CO*¥*2)) 
D{N})=-1.20 
E(N)= (H**3)*( BB-RE*F1I(N)*F120N)) 
64 CONTINUE 
CALL SGLU{A,BeCesD,ExsZZ) 
DG 67 N=2;3I1 1 
F2(N)= ZZ(N) 
67 CONTINUE 
F2(1)=0.0 
AFICKK)= F2{( 11) 
CF(KK)= BB 
ITF(DABS(F2CT1}3})¢ LTs 1-E-6) GO TO 72 
IF(KK.EGes1) GO TG 74 
IFC KKeEQ23) GO TC 23 
DXX= AF1(2)-AF1(1) 
DY Y=CF(2)-CF(1) 
BB=CF (1 )—AF1(1 )*DYY/DXX 
G6O TO 75 
74 §688=-88+.2005 
GO TO 75 
72 WRITE(6,73) 
73 FORMAT( 7014,9Xs°F2* 515Xs*F21%s15Xs*F22% 5 10Ks'F23") 
B21MM)=BB 
CALL DDETS (H»F2 »ZsNDM, IER) 
0G 78 N=1,11 
78 F21(NJ9=Z1N) 
CALL DDETS {(HsF21eZsNDM, IER) 
DG 79 N=1,51TI 
79 F22{N)I=Z(N) 
CALL DDETS {HeF223s3Ze2NDMZTIER) 
DO 80 N=1,iI 
80 F23(NJ=Z(N) 
DO 76 N=1,i1I 
WRITE(6054)F2(N) sF21(N) sF22(N) »F230N) 
76 CONTINUE 
CK RK CALCULATICN FOR EGe NOs 4 
READ(5531) B3{MM) 
BB=B3(MM) 
KK=0 
95 WRITE(6.,82) BB 
82 FORMAT(*0155Xse '8B3="s F106) 
KK=KK+1 
WRITE(6233) 
DO 84 N=3,II!1 
AUN)= (12-0-5*(H¥*¥2) *(RE*F OLIN) + (HART*#¥2)*(CO#¥*2))) 


B4 


87 


94 


92 
93 


98 


99 


100 


96 
C#KX 


ris 
102 


T00 


B(N)= ((CH*¥*3)*3*RE*FO2(N)—-3.-) 

CINI= €324005*(H**2) *¥( REXFOL(N) FTHART#*2) *¥(CO**#¥2))) 

DIN) =-1.0 

ECN )=(H*%*3 ) *(BB-RE#*(F1(N)*F22(N)42*F2(N)*#FI2(N)—F11I(N) 
*F21(N))) 

CONT INUVE 

CALL SOLU(AsBsCsDsEsZZ) 

DG 87 N=2,I1!1 

F3(N)=ZZ(N) 

CONTINUE 

AFI(KK)= F3CIT) 

CF(KK)= BB 

IF(DABS(CF3C(TI)) -« LT. 126-6) GO T0 92 

IF(KK.2EQ21) GO TG 94 

IF(KK.eEQe3) GO TO 23 

DXX= AF1{(2)—-AF1(@1) 

DYY=CF(2)-CF{1) 

BB=CF(1 )—-AFi(1)*DYY/DXxX 

GO T0 95 

6B=88+.005 

GO TO 95 

WRITE(63;93) 

FORMAT €* 0% s9Xs *F3* .15Xs47F3B14 s15X-.*FSB2° 5 10Xe4F33*) 

B3(MM)=8B 

CALL DDETS €CHeF3 »ZsNOMZI1ER) 

DG 98 N=1;3I!I 

F31(N)= ZN) 

CALL DDETS (HsF31+ZsNOM, ITER) 

DO 99 N=1>,5IT 

F32{(N)=Z(N) 

CALL DDETS (H»sF322Z+NDM,5 IER) 

DG 100 N=1,1I 

F33(N)J=ZI{N) 

00 96 N=1,1I 

WRITE(6,54)F31N) sF31(0N) sF32(N) sF33N) 

CONT INUVE 

CALCULATION FOR EQGQe NOs 5 

READ(5,31) B4{MM) 

BB=B4(MM) 

KK=0 

WRITE(6.102) BB 

FORMAT(* 0%s5Xs "84=",F1046) 

KK=KK+1 

WRITE(6s33) 

DG 104 N=3,II 

ACN)= €10-02e5*(H*¥*2) *( 2*REXF OLIN) +0 HARTE*2 )*(CO**X2) )) 

BIN)J= (4%*(H**3)*RE*F O20N)-340) 
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107 
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116 


135 
122 


CINI= (3440 .5%*(H*¥*2 )*( 2¥REXFO1L(N 94+ (CHART##2) *(CO*%2))) 


DIN)= —-120 


IOI 


E(N)= (H**2)%*( BR—-RE*{(F1(N)*F32(N) 4#2*F20N )*F22(N )4+3%*F 3{ 


x N})*FL20N)-—2 
Le¥FLI{N)*FS1L(UNI—-F21L(N)#F210N))) 
CONTINUE 
CALL SOLU{A+B8s,Cs,DsEsZZ) 
CO 107 N=2,I11 
F4(N)= ZZ(N) 
F4(1)=0.0 
AFI(KK)= F4C1I1) 
CF(KK)=BB 
IF (DABSCF4(1TI))- LTe 126-6) GO TO i112 
IFC KKeEGsi) GO TO 114 
IFC KK.2EQ23) GG TC 23 
DOXX= AF1(12)—-AF1(1) 
DY Y=CF(2)-CF{1) 
BB=CF(1)—-AF141 )*DYY/DXX 
GG TO i115 
8BB=-BB+2005 
GO TO 115 
WRITE(6,116) 
FORMAT( * 07 39Xs"%F 4% 515Xs *F41%9515Xe*F42%*5 10Ks"*F43% 5 
B4(MM)= 8B 
CALL DDETS (HsF4 +ZsNDMSZIER) 
DG 117 N=1,i1 
F41(N)= ZN) 
CALL DDETS (H»F412Z+NDM, IER) 
DO 118 N=1,ilI 
F42(N)= ZN) 
CALL DDETS (HsF42eZ+NDM, TER) 
DO 119 N=1,II 
F43(N)= Z{N) 
pO 120 N=i,II 
WRITE(6,54) F4(N),F414N) sF42(N) »F430N) 
CALCULATICGN FOR EGe NOs 6 
READ(5531) BS(MM) 
BB=B5(MM) 
KK=0 
WRITE(6,122) BB 
FORMAT(°01,5X» "85=",F1026) 
KK=KK+1 
WRITE(64,33) 
DG 124 N=3,11 
A(N)=(1-065*(H**2) *( B*¥RE*FO01(N)4(HART*®*2)*(CO**2))) 
BIN) =(CH**3) XS *RE*XFO2(NI)—3e ) 
CON) =(30405*(H*¥*¥2 )*¥(3*REXFO1(N) 4+ (CHART ¥*2 )*(CO**2))) 
DIN)=-1.0 
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C#** 
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EQ(N)= (H**¥3)4(B8- RE*( FIN )*F421N)42%F2(N) *F32(N)43*F3 


(N)*F220N)+4 


1 4*F4(N)4F12(NI-—3¥F11(N) *F410N)—3*F210N)*F31(0N))) 
124 CONTINUE 


* 


CALL SGLU(AsBsCsDsEsZZ)} 
DO 127 N=2,11 

FSCN)= ZZN) 

FS5(1)=0.0 

AFICKK)= FSCI1) 
CF(KK)=BB 

TF(DABSCFS( TID). tT. 12E-6) GO TO 132 
IFC KKeEQe1) GO TG 134 
IF(KK.EQe3) GO TG 23 
DXX= AF1{2)-AF14@1) 
DYY=CF(2)-CF(1) 

BB=CF (1 )—-AF1(1)*DYY/DXX 
GG TO 235 

BB=BB+.005 

GG TQ 235 

WRITE(62136) 


FORMAT ('0* 59Xs4F5S*%215Xs*FS1%s15Xs*F52%5 


B5(MM)=BB 

CALL DDETS (HsF5S »ZsNDMsIER) 
DO 137 N=1,I1I 

F5S1(N)=Z(N) 

CALL DDETS {(HeFS1lsZsNDM, IER) 
DO 138 N=i1,I11 

FS2(N)= ZiN) 

CALL DDETS (HsF520ZsNDMs IER) 
DO 139 N=1,11 

FS3(N)= ZEN) 

DO 140 N=1,11 

WRITE(6:54) FS(N)seFS1(N) »FS2(N) sFSSIN) 
CALCULATION FOR EQNse NOs 7 
READ(5.231) BE{ MM) 

BB=B6(MM) 

KK=0 

WRITE(6s202) 8B 
FORMAT(*0%55Xs *B6=% 5F 1026) 
KK=KK+1 

WRITE(6,33) 

DO 204 N=3sII1 


10Xs*F53*) 


AGN9=(1 20-0 0 S¥(H¥*2 *( 4*REXFOL(ND + CHART R*2) X(CO*KK2) ) ) 


BIN)J= (6 *(H**3) #RE*#FO2(N)I-3-0) 


C(IN)= (340 e5#CH**2 )*(4*RE*XFOL(N)I FM HART ## 2) *&(COX*2) ) ) 


DIN)=—-1. 


E(N)= (H*¥*3) *(BB-RE*{F1(N) *FS20N)4+2*F20N) *F42(N)43%F3 ( 


N)*F32(N) 


204 


207 


214 


4i5 


97 


)*F41(N)—-2* 


2F31(N)*F311N))) 


CONTINUE 

CALL SOLUTAsBsCs0sEsZZ) 
DO 207 N=2,I1I 
F6OUINJ=ZZUN) 

F6(1)=0-0 
AFICKK)=FE6CII) 
CF(KK)=88 
TFC(DABS(F6(1I1I))« LT. 12E&-6) SO TO. 2t5 
IFC KKeEQe1) GO TO 214 
IF(KK2EQ.e3) GO TG 23 
DXX= AF1i(2)-AF1(1) 
DYY=CF(2)-CFi1) 
BB=CF(1)—-AF1(1)*DYY/DXX 
60° TO 201 

2B=BB+.005 

GG TO 201 

WRITE(6,216) 


FORMAT (* 0" 59Xe*F6% 515X%s "F661" s15Xs "F622" s 


BS(MM)=BB 

CALL DDETS (HeF6 «ZsNDMZIER) 
DO 218 N=1,II 

F6ELIIN)J=Z0N) 

CALL DDETS {HesF61sZ2NOM,s IER) 
DG 219 N=1,I1I 

F62(NJ=ZIN) 

CALL DDETS {(H»F620ZsNDM, ITER) 
DO 220 N=1,5I11 
FO3Z(0N)J=Z(N) 

DO 222 N=1,II 

WRITE(6254) F6E(N)sFE1(N) »-F62TN) »FOS CN) 


CALCULATION FOR PRESSURE DISTRIBUTICN 


WRITE(6:97) BO(MM),B14MM)5B82(MM) .B3(MM) 


»B6O(MM) 


FORMAT(*° 0% ,5X%s*BO0=*5F1026/ "Bl=*" 5F10e6/ 


B3=*+F10.6/ 


CUR=0-0 
DO 800 J=1,21 
DEL(J)I= (120-02025*(J-1)) 
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1+4*F4(N) *F22(0N) 45 *F SN ) *¥F1L20N)—-4*F 1 LOND *¥F510N)-4*F210N 


10Xe*F63" ) 


»B4(MM),85(MM) 


*"B2="*,FL00«6/ *¢ 


1SXs "84=",F1046/ *85=*,F10-6/"*B6='.F1026) 


EE(I)= (44/30) #HART*¥ (COZ (DEL (J) )4C14C2*(DEL( J) )4+C3*((D 


EL(J))**2)) 


*2))% 


2 F0(41)*(CO0*1.38E3)) 


1*(CUR+2*HART*(CO/S( DEL (J) )4+C14C2*( DELI) )4+C3*( (DEL( J) )* 


— i 


* 


( 


{ 
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PGRAD(JI=1-—2e/(REX(IDEL(I)) **3))* (800MM)4B1(MM)4DEL7 J) 

* +#B2(MM)* 

1( (DEL J) )**2)4B3(MM)*( (DEL (J) )**3) 4B4(MM)*( (DELS) ) ** 

%e 4)+ BS(MM) * 

2(¢(DELCIJ))**5) +B6(MM)#( (DEL (J) ) **6) I-TEECII/SRE)D 
800 CONTINUE 

S=0205 

ND=21 

CALL DGSF(S»PGRAD+sZsND) 

WRITE(6s810) 
810 FORMAT(*0%512Xs *DEL*s20Xs*PGRAD*»s20Xs*PRESS!* ) 

DO 840 J=1i,21 

WRITE(63820) DEL{J)sPGRAD(J).Z{(JI) 
820 FORMAT(470*,5XsF10206915XsF10063915Xs9F1026) 
840 CONTINUE 

CALL DGSF{(S.ZseZPseND) 

WRITE(6s830) ZPT21) 
830 FORMAT (40°,5X%e *LOAD="*,F10.6}) 
23 CONTINUE 

STOP 

END 


SUBROUTINE FCT {XsYsDERY) 

DOUBLE PRECISICGN Xe ¥(3) sDERY(3) »CB,FO0(55).F01(55),F 02 
5 (55), 

1 FO3{1(55) sREsHART».CB2TO 

COMMON FO>sFO01eF023F03%RE 2sHART2CBsCO0 IT 

DERY(1)=Y¥(2) 

DERY(2)=Y (3) 

DERY(3)= CB-Y{(2) *Y(2) *RE+( HART**2) *(CO**2)*Y( 2) 

RETURN 

END 


SUBROUTINE GUTP(XsYsTCERYsIHLFesNDIMsPRMT ) 
DOUBLE PRECISION XeY(4)2DERY(3) 0 PRMT(5) sREe®HART »CB2CO 
DOUBLE PRECISION ETAINT(55 )sF0(55)5F01(55) »sF02(55) »FO 


*ENSC5SS) 

COMMON FOsFO1sFO02sFO03sRE sHART»CBsCO2IC 
Ic=1C+1 

ETAINT (1C) =X 


FOC(IC)=yY(1) 
FOI(C IC)d=Y(2) 
FO2(IC)=¥(3) 
FO3(1TCI= DERY(3) 


{ ) a+ ) ) eC TE ee (OL) i I\eS~—}={L) é a : 
bas : A i Me is ae or puta 2 : > 
i Ge } (MM } + ft ((G) 530) >) ( pEHe (SKE ECL DL rr! 
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WRITE(6,23) ICsETAINTCIC) sFOCIC) sFOICIC) sFO2(IC) sFO3(I 
Cc) 

FORMAT( *0*,5X%s12s2Xs5F1526) 

RETURN 

END 


SUBROUTINE SGLU(AsB»sC2D2EsZZ) 
DOUBLE PRECISICN A(55)38(55),01(55).D(55),E(55)4.2Z2(55) 
»BOT( 55) 


1TOP(55)5R(55),.CC(55) »ANUMEs ADEN 


COMMON /FIRST/ZIL 

TI=41 

IlI=11i-1 

IUPPER=II-2 

R{(2)=0.225 

CC(1)=00 

IFCTIL*EQe21) GO TD 20 

CC(2)=0.20 

GO TG 30 

CC{2)= —-02025/4.0 

DO 300 N=3,ii!I 

BOTIN)= 81{N)4C(N)*RIN—1)4D1N)9*%RIN—-1 9 *RIN—-2) 

R(N)=-A{N)/BOTI(N) 

TOP(N)= —-(CIN) *CCIN—1)4D00N)*RON—-2) * CCI N—-1) 4+D0N) *CC(N-2 
))4+E CN) 

CC(N)= TOPIN)/SBOT(N) 

CONTINUE 

WRITE(6.36) 

FORMAT (90° ,12Xe"RN%215Xs *CCN® ) 

ANUME=CC(11-1)-CC( T1i-2)7(42-R( 11-2) ) 

ADEN=3e/(4e-R( 11-2) )-R( 11-1) 

ZZ(LI1)=ANUMESADEN 

DOD 39 N=1,I1UPPER 

I=TI-N 

ZZ( 1P=RCTI*ZZ(14+1)94CC(T) 

CONTINUE 

RETURN 

END 


cy 
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* IN THE SIXTH CGLUMN INDICATES 
CONTINUATICN FROM PREVIOUS LINE 


HK RRR KKK AR RRR KR KH EK 


C##** FORTRAN PROGRAM FOR CHAPTER 4 
C#*** THE PROGRAM CALCULATES PRESSURE DISTRIBUT IONsL OAD 


Cc CAPACITY .sVELOCITY DISTRIBUTIGN»sFRICTION FACTOR FOR 
c MHD SLIDER BEARING 
CxX#e*xe SUBROUTINES FROM IBM—SYSTEM/360 SCIENTIFIC SUBROUTINE 
c PACKAGE ARE USED IN THIS PROGRAM 
C#e#K* THE APPLIED MAGNETIC FIELD MAY BE UNIFORM OR 
aS LINEARLY INCREASING TYPE 
DOUBLE PRECISION HART(5) »FLOW(5).FLO15),SY8(21) sHBAR( 
* P2l)s 


IPGRAD{( 21) sPRES(2)2Z121)3ZP{21) »~DELXs,DELYsDFLOW.»A121}).8 
* aCaDoPs 
ZALPHATL213,B8Y »RAT(7}).,DT(7) »sBB»,D0D(5) »sPLi2ijeDYs,Cilel, 
ak C2{21), 
3B €3(21)0UL 21221) DUY(21221),5Y(21.21) DTO(21)227T(21) 
C#**K* FOR UNIFORM MAGNETIC FIELD 111.GT 01 
Cx#x##e* FOR LINEARLY INCREASING MAGNETIC FIELD LLi-i 


Lii=2 
C*¥*** READ HARTMANN NUMBER, MASS FLOW AND FILM THICKNESS 
c RATIO AT INLET OF BEARING 


DO 10 L=1-3 
READ(5.11) HART(L).FLOIL) »DDO CL) 
11 FORMAT(3F1026) 
WRITE(6s12) HART(L) 
12 FORMAT(* 0" 45X~s "HARTMANN NUMBER=" »F1026) 
FLOW(1)=FLO(L) 
Cx#kK* MESH SIZE 
C#¥#x*%* CALCULATE HBAR ALONG THE LENGTH GF BEARING 
DELX=.05 
DELY=.05 
DFLOW=0202 
DO 31 IK=1,4 
RATCTIK)=DD(LI+02010*(IK-1) 
DTCIKI= (RATCIKI-120)/2020 
WRITE(6332) RATCIK),O0T(IK) 


a2 


CHEEK 
CRE X 


100 


iS 


CKKEX 


20 


FORMAT(°04sSXs *"HBAR AT INLET="*sF10e62"0%s5Xs *OT=",F il 
% 024) 

CALCULATION ARE CARRIED CUT USING THE ASSUMED MASS 
FLOWsLATER REPEATED WITH THE IMPROVED MASS FLOW 

KK=0 

KK>KK41 

DO 33 I1=152 

WRITE(6.15) FLOW{I1) 

FORMAT(*0*,5Xs *MASS FLOW=*,F10.6) 

iF (LLL.EQ.1 )} GG TO 101 

CALCULATION IF THE APPLIED MAGNETIC FIELD IS UNIFORM 
DG 20 1=1,21 

HBARC(T)I= (CRATCIK )-DT( 1K )¥(1-1)) 

ALPHA I) =( CHBARC I) )**2)*( (CHART (CL) )**2) 

SYBCT)= FLOW(TID/SHBAR(I) 

P=( ALPHA( I) )**025 

PL(I)= (ALPHA I) )**0.5 

B= (DEXP{—P)4DEXP(P)-22-0) 

C= (-SYB(1)*P*( DCEXPI—-PJ—DEXP(P))/B -1-0) 

D= 120-DEXP{(—P )+{(DEXP(—P)-DEXP(P) )¥*112—-DEXP(P) 4P#DEXP{ 
* PJISIDEXP(P) 

1*8) 

ALI) =(C/D)*ALPHATI) 

C3C1I)= (DEXP(PIZ( P¥{TDEXP{—-P)-DEXP{P)))) *(1-4+(AT TI) 70P 
*% *ke2)9-(ATTIS 

1¢(P**2)*DEXP(P)))) 

C2CT)9= (12/0 P#DEXP(P) ))*( C301) *P*¥DEXP(-P)4+( ATT) /(P**2) 
ae 2 

Creer -c2eri—-catti 

BB= (HBAR(1)41.0) 

PGRAD(1)=12-/8B) *HBAR(I)* (CTHART(L) )**2)*SYB(I ) +A) 
* / {(THBAR(T)) 

1*x*x2) 

CONTINUE 

GO TO 102 

CALCULATION YF THE APPLIED MAGNETIC FIELD IS LINEARLY 
INCRASING 

FOR LINEARLY INCREASING CASEs,THE MAGNETIC FIELD AT 
OUTLET IS Twice THAT OF FNLET 

DO 30 {=1,21 

BY=(1.+.-05*{(I-1)) 

HBAR(I)= (RATT IKI-DTCIK)*(I-1)) 

ALPHA (1)=( CHBARCTI ) ) #42) *(THART(L ) ) **2)*( BY*X2) 

SYB(1)= FLOWC(ITID/SHBARCT) 

P=(ALPHA( I) )**025 

B= (DEXP(—P)4DEXPIP)—-220) 

PL(I)= (ALPHATI))**0.5 

C= (—-SYB(1)*P*#*( DEXP(-—P)-DEXP(P))/8B -1-0) 
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30 
CK 
102 


490 


43 


BS 


CK K 


39 


42 
41 
31 
CHEK * 


* 
1 


= 


* 


oe 
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D= 140-DEXP{(—P )4+{DEXP{(—P)—-DEXPCP) 2*(1.-DEXP(P )4P*DE XP ( 
PIISTIDEXP(P) 


*B) 


ACT) =(C/D) *ALPHA(I) 


C3(1I)= 


(DEXPTPI/¢ P*(DEXPI—P)J—-—DEXP(P)))) *€1.4¢ACTT) 70P 


*¥*2)9-C ACIDS 
1((P**2)*DEXP(P)))) 


Cc2(1j= 
2) 


C2-/(P#¥DEXP(P)))*(C3( 1) *P¥DE XP(-—P)4( AC TIS (P#*2) 


Citib= —-Ceci)-c3 41) 

BB= {HBAR{1)+1.0) 

PGRADCI)= (26/88) *HBARC 1) *(CHARTC(L) ) **2)*SYB( 1) *B8Y*1.45 
+ATTI/ (CC HBAR 

1€1))**2) 

CONT INVE 

CALCU FOR PRESSURE DISTRIBUTION AND COAD CAPCITY 


ND =21 


CALL DGSF(DELXsPGRADsZ+ND} 
PRES(11)=Z(21) 

WRITE(6340) 

FORMAT(70%.s10Xs *PGRAD * »18Xs*PRESS'*) 
CALL DCOSF (DELXsZsZPsND) 

WRITE(6243) ZP{21) 

FORMAT(*0%,5X»s "LOAD=" 3sF10-6) 

IF {DABS(2Z(21))eLT2»-0005) GO TO 39 
FLOW(1T1+1)= FLOW(IT)+DFLOW 

CONTINUE 

IF (KK.2EQ24) GO TO 39 

BETTER APPROXIMATION ON FLOW 

DX=PRES (2)-PRES(1) 

DY=FLOW(2)-FLOW(1) 
YL=FLOW(1)-(PREST1)*DY/DX)} 

FLOW(1)= YL 

GO TO 100 

DO 41 1=1,21 

WRITE(6242) PGRAD(I).Z(1) sHBAR(I) 
FORMAT €20445XsF10e6015XsF1006310XsF 1046) 


CONTINUE 

CONTINUE 

CALCULATE STREAM FUNCTION DISTRIBUTION 
NP=21 

DO 152 1=1,21 


SY(1Ts1)=0.20 

SY(1s21)=SYBC1) 

DO 154 J=2,20 

DY=DELY*(J-1) 

SY(1IsJ)= C1I(1)4C201V*¥DEXPCPL II )*OY) 4+C3( T)*#DEXPC—PLII)* 


DY) 


~CACI)* 
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IDYISICTPLCI) )**2) 
154 CONTINUE 
152 CONTINUE 
WRITE(6,157) 
157 FORMAT(/S/5Xe"*STREAM FUNCTION DISTRIBUTION IN THE BSEARI 
ak NG *} 
DQ 160 I[=1,21 
WRITE(6:167) (SY@1I.J)sJ=15,21) 
167 FORMAT(*® 04 ,.2Xs21F623) 
160 CONTINVE 
C#¥*x** CALCULATE VELOCITY DISTRIBUTION IN THE BEARING 
N=20 
C#*** CALCULATE *U* AND *DUY*® AT EACH SECTION OF THE BEARING 
DG 170 f=1,3,21 
UL Is,ij)=1.0 
ULIs2)= (SY(1,5)-64*S5SY(154)418*SY(1Is3)-10*SY(13,2)-3*SY( 
* Ts1)9/7(12% 
iDELY) 
UtLIsNi= (SYCIsN-3)-6*SY( IsN-2)418*SY(C1I5N—-1)—-10*%SY(14N) 
* -—3*eSY(I,N41 
1))/{1-12*DELY) 
U{I.2N4+13=0.0 
DUY(Is2)= (112*S5Y(155)-S6*SY(1,4)4114*%*SY(153)-104*SY(I, 
% 2})+35*SY(1>» 
1 #1))/7¢€12*(DELY**2)) 
DUYT1Ts2)=(-SYC1,5)44*SY(154)46*SY(153)-20*SY(152)411*S 
% Y(Iisi1d)/ 
I{L2*{(DELY**2)) 
DUY@I»sN)I= (11%*SYCI19N4+1)—-20*5Y(C19N) ¢6%*SY(1,N-1344%*SY(15 
* N—-2)-SY( 1s 
IN—3))/012*{ DELY**#2)) 
DUY (12N41)=( 35*SY(1.4N41)-104*SY(15N)4114*SY(1I.N—1)-S6* 
% SY{IsN-2)+ 
1 LI*SYTIsN-3))/012*(DELY**#2) ) 
DO 172 J=3219 
U(IsJ}= (SYL1s J-2)-8*SYC1,I5-1)48*SY(1554+1)-SY(1s542))/7 
* (12*DELY) 
DUYCIsJ)= (-SY01434#239416%SY(1.J54+1)—-30*SY( 1,45) 416*S5SY(I. 
x J-13-SY(1I» 
1 J-2))/(12*(DELY**2)) 
172 CONTINUE 
170 CONTINUE 
WRITE(60177) 
177 FORMATLS/5Xe"U DISTRBUTION IN THE BEARING! ) 
DO 178 {=1,21 
WRITE(60179) (UCIsJ) sJ=1521) 
179 FORMAT(* 08 »2X%e21F623) 
178 CONTINUE 
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WRITE(6,180) 

180 FORMAT{I/Z/5X-*DUY DISTRIBUTION IN THE BEARING® ) 
DO 182 11:21 
WRITE(62185) {DUY(I. J) sJ=1321) 

185 FORMAT(*0',2X%s21F6.3) 

182 CONTINUE 
DG 200 I=1>,2i 
DTOCI) =DUYC1,1)4*(1-4¢/HBARII)) 

200 CONTINUE 

C##*e* CALCULATE SHEAR FORCE 

DOB 202 {=1.21 
WRITE(6,203) OTOC(I) 

203 FORMAT(*04,5XeF E623) 

202 CONTINUE 
CALL OGSF(DELXsDTCsZIsNP) 
WRITE(6.,201) ZIiI({NP) 

201 FORMAT{(*04,5X. "SHEAR FORCE =*,F623) 

10 CONTINUE 

STGP 
END 
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* IN THE SIXTH COLUMN INDICATES 
CONTINUATICN FROM PREVIOUS LINE 


Se RK RR RK gO RK RK KK KK K 


Cx#x*** FORTRAN PROGRAM FOR CHAPTER 4 
Cx*k*e THE PROGRAM CALCULATES PRESSURE DISTRIBUTIONsLOGAD 


c CAPACITY s VELOCITY DISTRIBUTIONsFRICTION FACTOR FOR 

€ MHD SLIDER BEARING 

C**** SUBROUTINES FROM IB8M-SYSTEM/360 SCIENTIFIC SUBROUTINE 
c PACKAGE ARE USED IN THIS PROGRAM 


C#¥#%kxe THE APPLIED MAGNETIC FIELD IS STEP TYPE 
DOUBLE PRECISION HART(5)sFLOW(5)eFLO(S5).SYB(41)sHBAR{ 


% 41}> 
IPGRAD(41 ) sPRES(2)0Z141)4ZP(41) sDELXK»DELYs, Al41) »eBoCsDs 
x Pe 


ZALPHA(41}),8Y sBBeBETA,DFLOW( 5) ».PGR(41)52Z2(41),RAT(7).D0 
* T(7},.DB 5) 
3 DHS) sPL(41)2C1041).C2741) «C3(41) sDVsSY(141 521 )5U(4142 
x 1) sDUY{(41 5,21 
4),0TO(41),.,Z1(41) 
C*¥*#*x* READ HARTMANN NUMBERSMASS FLOW » INCREANT IN MASS FLOW 
Cx**%* STEP LOCATION AND FILM THICKNESS RATIO AT THE INLET 
DO 10 J=ie1 
READ(5,11) HARTCJ)sFLOII) sDFLOW( J) 20815) »DH (JI) 
11 FORMAT{(5F1027) 
WRITE{6212) HART{(J) 
12 FORMAT( 407 .5Xo "HARTMANN NUMBER=* .F1026) 
FLOW(1)=FLO( J) 
NP=41 
CkxkKE* MESH SIZE 
DELX=02025 
DELY=.05 
Cekekx*x BETA DIMENSION IS STEP LOCATION GF APPLIED MAGNECTIC 
i FifTwvlo 
DO 21 IJ=1,53 
BETA=D8(J)+02025*(1J-1) 
NI=BETA*4025 
NN=NI41 
WRITE(6317) BETA»NN 


17 FORMAT(°0*,5Xs *BETA=*4F1002695Xs *NN="512) 


C#x*e*e CALCULATE HBAR ALONG THE LENGTH OF BEARING 


DO 31 IK=isi1 
RAT CIK J=DH(C J)—-0205%*{ IK-1) 
DTCIKI= (DABSTRATCIK)—-120))/74020 
WRITE(6.s322) RAT{L IK) »sDT(IK) 
32 FORMAT(*0*,5Xs*HBAR AT INLET=" .F10260°0%s5Xs*DT='"5 
x 26) 


C#¥*** CALCULATION ARE CARRIED CUT USING THE ASSUMED MASS 


S 


FLOWsLATER REPEATED WITH THE IMPROVED MASS FLOW 
KK=0 
100 KK=KK+1 
DO 33 11=1,2 
WRITE{6,15) FLOW(II1) 
15 FORMAT{(*0%.,5Xe *MASS FLOW="5F1026) 


CHEEK CALCULATION FOR NO MAGNETIC FIELD PART 


DO 13 N=1.NN 

HBARIN)I= {RATCIEKI-DTCIK) *(N—-1)) 

SYB(N)= FLOWCIETISHBAR(N) 

PGR (N)= (66-12e*SYBIN) 9 /(THBARIN) )**2) 
13 CONTINUE 

BB=0 20 

ND=NN 

CALL DGSF( DEL XsPGRsZ»ND} 


C¥x** CALCULATION FOR MAGNETIC FIELD PART 


L=0 

DO 20 I=NN»NP 

HBAR{I)V= (RATT IKI-DTCIK})*(1I-1)) 

ALPHA (1 )=( (HBAR(1))* 42) *{(THART( J) )**2) 
SYBCT)= FLOWCITI/SHBARTI) 

P=( ALPHA( 1) )*#025 

PL(TJ= (ALPHATTI))**025 

B= (DEXP(—P)4+DEXPIP)-2-0) 

C= (—-SYB(1)*P*( DEXPI-—P)-DEXP(P))/7B —-1.0) 


112 


FiO 


D= 120-DEXP(—P )4+(DEXP(-—P)-—DEXP(P) 9*(1.2-DEXP(P)4+P*DE XP( 


* P)I/S(IDEXP(P) 
1*B) 
ACI) =(C/DI*ALPHAII)D 


C3(1)= (DEXP(P)/( P*X(CDEXP(—P)-DEXP(P)))) *(1-4+L ACT) /(P 


* &*#2)9-CATI)S 
1((P**2)*DEXP(P)))) 


C2C1TI= (1-6/1 P¥*DEXP(P)))*(C3( 1) ¥P*¥DEXP(—P)4( ATT) /(P**2) 


* )) 
CLO T}= —C2¢1)-C301) 


PGR (1)=(2.*8B) #HBARCI)*® (THART( J) )**2)*SYBCT ) +ACT) 


* / (THBARIT)) 
1**2) 
L=L+l1l 
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PGRAD{LI=PGR{I ) 

WRITE(6245) A@l1),C3(1),C2(1)4.C1(1) 
FORMAT (404% ,5X,4F10.6) 

CONTINUE 

ND=NP41—NN 

CALL DGSF (DELXsPGRADsZZ7»ND?} 
NM=>NP—-AN 

DO 25 IP=14MM 

Z(NN4iIP)= Z(NN) 4+ZZ{(1P41) 
CONTINUE 

PRES( 11 )=Z( NP) 

ND=NP 

CALL DGSF {DELXsZs2PsND) 

IF (DABSTZINP)).LT.~200005)G0 TO 39 
FLOWCII41)= FLOW( TID 4+DFLOW{ 53> 
CONT INUVE 

IF (KK.EQ.3) GO TO 39 

BETTER APPROXIMATICN CN FLOW 
DOX=PRES(2)—PRES{1) 
DY=FLOW{2)—-FLOW (1) 
YL=FLOW{(1)-—fPRES(1)*DY/DX)3 
FLOW{1)= YL 

GD TO 100 

WRITE(6.,37) HBAR{INP)>,ZINP) 
FORMAT €'0% 5X sF102026%010X3,F10 26} 
WRITE(6543) ZPUNF) 
FORMAT(70*,5Xs "LOAD=" »F1026) 
CONTINUE 

CONT INVE 

CONTINUE 

CALCULATE STREAM FUNCTION ODOFTSTRIBUTIGON 


DO SO 1=1,NN 

SY(1,1)=0-20 

SY€1.21)=SYBUT) 

DO 52 J=2220 

DY=DELY*(J—-1) 

SY(1IsJ3)=(12-2*SYB( 1) D*DY#%3 -—-(2-.-3*SYB( 1) )*DY**2 +DY 
CONTINVE 

CONTINUE 

DO 152 I=NNo41 

SY¥(1.1)=0-0 

SY¥(1.21)=SY81) 

DO 154 J=2.20 

DY=DELY*{J-1) 

SY(T,J)= C1(194+C27T1I)*DEXPI(PL(I1)*0Y )4C3(1T D*¥DEXP{—PL(1)* 
* DY) —-CACI)* 

IDYIS (CPL) )**2) 
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154 CONTINUE 
152 CONT INVE 
WRITE(6,157) 
157 FORMAT(IS/5X.*STREAM FUNCTION DISTRIBUTION IN THE BEARI 
* NG *) 
DG 160 1=1,41.2 
WRITE(64167) (SY(1IsJ)sJ=1,21) 
167 FORMAT{(70%,;2X.,21F623) 
160 CONTINUE 
CHEK CALCULATE VELOCITY DISTRIBUTIGN IN THE BEARING 
N=20 
C##eX* CALCULATE "U*® AND "®DUY* AT EACH SECTION CF THE BEARING 
DO 170 1=1.41 
U{1I,1)=1.0 
ULI.s2)= (SY(1,.5)-6*SY(124)418*SY(1.3)-10*SY(1.2)-3%*SY{ 
x 11) )/(12* 
1DELY) 
U{I~s~N= (SY(IsN—-3)—-6*SY(15N—-2)418*SY(1IsN—-1)-10*SY(1I.N) 
* —~3*SY(IsN4+1 
1))/7{-12*DELY) 
UCI»sNt1)=0.0 
DUY(CT21)= (11*SYC1,5)—-S56*57(154)4114*S5SY(1s3)-104*SY( Is 
am 2)4+35*SY( 1» 
1 1))/7(12*(DELY**2)) 
DUY(1,22)=€-SY(125)44*SY(01,4)+6*SY(123)-20*S5SY(1I52)411%S 
* Y(I21))7 
1(12*(DELY**2)) 
DUY(CTsNS= (11%*SY(15N419-204SY(15N)46*SYCI2N-1)44*SYCUIs 
a N-2)-SY(I; 
IN-3)9)9/(12*(DELY*%#2)) 
DUY (I1sN+1)={ 35*SY(C TsN41)—-104*SY(1,N)4114*SY(1,N—-1 9-56*% 
* SY{IsN-2)+ 
1 11*SYCIeN-3))/7(12*(DELY**2)) 
DO 172 J=3219 
UlIietJ)= (SYC1,I5-2)-8*SY( 1s 3-1) 4+8*S5SY(1,541)-SY(1sJ4+2))97 
%K {12*DELY>) 
DUYCI,J)= (-SYC1.542)416%S5SY(1554+1)-S30*SY( I,J F16*SY(I5 
* J-1)-SY(i>+ 
1 J-2))/012*(DELY*#*2) ) 
172 CONTINUE 
170 CONTINUE 
WRITE(6s177) 
177 FORMAT(//5X%e1U DISTRBUTION IN THE BEARING*) 
DO 178 1=1%41s2 
WRITE(64179) (UC TsJI) eJ=1521) 
179 FORMAT( *0%52Xs21F623) 
178 CONTINUE 
WRITE(6.,180) 
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185 
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CK 


200 


203 
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201 


FORMAT(//5Xs"*DUY DISTRIBUTION IN THE BEARING® ) 
DO 182 I1=12413,2 

WRITE(62185) (DUYTIs, JJ eJ=1521) 
FORMAT(*0*,2X,21F623) 

CONTINUE 

CALCULATE SHEAR FORCE 

DOGO 200 1=1,41 

DTO(i) =DUYC 1,1)9%*(1-/HBARCI)) 

CONT INUE 

DO 202 1=1,4i42 

WRITE{(6.203) DTG(I) 
FORMAT(* 0% 5 SXeF 623) 

CONTINUE 

CALL DOSF(DELXsDTO.Z15NP) 
WRITET6,201) ZI{NP) 

FORMAT (°0*,5Xs *SHEAR FORCE =*.F6e3) 
STOP 

END 


~ EE 


CHKK* 


CK IK ke 


LO 


FH HC ie eR Fe 3 ak I eK ak ak a ee a aK a oe ok 


* IN THE SIXTH COLUMN INDICATES 
CONTINUATION FROM PREVIOUS LINE 


BE EK He oe ea ae ek a he fea eae a a BK a OK a aK aK 3 ee 


FORTRAN PROGRAM FOR CHAPTER 5 

THIS PROGRAM CALCULATES POTENTIAL DISTRIBUTION FOR 
FINITE WIDTH PARALLEL PLATE BEARINGsFROM POTENTIAL 
DISTRIBUTIGN ELECRRIC FIELDsCURRENT DENSITY,CURRENT 
DENSITY GRADIENTS ARE CALCULATED 

THE POTENTIAL EQUATICN IN FINTE DIFFERENCE FORM IS 
SOLVED USING POINT SUCCESSIVE-OVERRELAXATION METHOD 
SUBROUTINES FRCM IBM-SYSTEM/360 SCIENTIFIC SUBROUTINE 
PACKAGE ARE USED IN THIS PROGRAM 

DOUBLE PRECISION AlsA2sA3sA42A52DXaDYs SYY*SYNsDIF sDI > 


* OMESERR, 


ISY€1014101)+.EX(101,101)5EZ(1015101) 0X1 sDV1 sA11,A225A3 


* 32A443A559 
2 $(101).S4N0101)5CUL{101)52Z0101) »CX(101%101).C2Z(1015,10 
* 1) 


DATA FOR CALCULATION 
M=65 

N=41 

M1i=M-1 

N1=N-1 

DX=0 20125 

ODY=0 20125 

OME=1.7 

IT=350 

ERR=0.200003 

LL=33 

LP=LL+1 

LK=49 

LO=LK+1 

LKK=LK-1 

LM=LK+1 

NL=21 

NLL=NL~-1 

NP=NL-1 

SY SYMBOL INDICATES POTENTIAL IN THIS PROGRAM 
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STARTING VALUES GF SY AT GRID POINTS 

DO 10 1=1;31K 

SYCI»N)=0-20 

DB 20 J=isNl 

SY(1IsJ)=0.23 

CONTINUE 

CONTINUE 

DDB 13 I=LO.M 

SY( 1I,NL)=0.20 

DO 14 J=14NLL 

SY(1TsJ)=0e1 

CONTINUE 

CONTINUE 

CALCULATE COGEFFICIENTS OF FINTIE DIFFERENCE FORM OF 
POTENTIAL EQUATION FOR SMALi STEP SIZE 
A1l=1./1DX**2) 

A2=1e/(DX**2) 

A3Z=1./(DY**2) 

A4=12/ (DY**2) 

AS5S= ((DX*DY)**2)/(24(0X**24DY* *2) ) 

DX1=0.025 

DY1=0.025 

CALCULATE COEFFICIENTS OF FINTIE DIFFERENCE FORM OF 
POTENTIAL EQUATION FOR LARGE STEP SIZE 
A11=1./(DX1*%2 ) 

A22=12/(DX1**2 ) 

AZB3B=1e/ (DY1%**2) 

A44=12/( DVI **2 ) 

ADS=((COX1I*OY1)**2)7 (2*(DX1**24+DVY1**2)) 

KK=0 

DI=0.0 

STARTING FROW INLET CF BEARING AND MARCHING TOWARDS 
RIGHT SIDE 

CALCULATICN CONTINUES UP TO THE END OF SIDE ELECTRODE 
USING SMALL STEP SIZE 

DG 45 1=1.LL 

SY(1.1)=2202 

Ti=1+1 

Ii1i=i-1 

DD 47 J=25N1 

JJ=J+1 

JJIJ=J—1 

IF(1-2EQe1) GO TO 48 

SYY= (AL¥*SY( ITs J) FA2*SY( TIT s JI FA3*SYC 1s JID +A4*SY( 15 JJ 
* ))*AS 

GO TO 49 

SYY= (AL¥*¥SY(CTIsJ)4A24*SYCII sJ)tA3*SVY( Ts JID FA4*SY( Is JIS 
* )*A5 
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SYN=OME*(SYY-SY( 125) )4+SY(12J) 

DIF=DABS(SYN-SY(13;J))/SYN 

DI=DMAX1 (DI.,OIF) 

SYCTsJ)=SYN 

CONT INUE 

CONTINUE 

CALCULATION FROM END OF ELECTRODE ONWARDS FGR 

INSULATED SIDE PLATE USING SMALL STEP SIZE 

DO 35 1I=LP.s.LkkK 

I{l=1+1 

Ili=iI-1 

DD 40 J=15Ni 

JJ=J+1 

JJJ=J5-1 

IF(J.-GT.1) GO TO 33 

SYY=(AI*SY( IITs JI F4A2¥*SY(C ILI s JF AS*SY( 1s 5S SI4A4*SY(T555)) 
*x AS 

GO TO 34 

SYY= (AL¥SY(C LI. jJ)4A2*SY( ILI s JV 4A3*SY( 1s JJ)FA4¥*SY( 13 JIS 
3) *AS 

SYN=OME*(SYY-SY(1,J))4+SY(14J) 

DIF=DABS(SYN-SY(IsJ))/SYN 

DI=DMAX1 (DIsDIF) 

SY(1I,J)=SYN 

CONTINUE 

CONTINUE 

CALCULATION AT TRE JUNCTION LK GF SMALL AND LARGE STE 
POSPZE 

DG 80 J=1lsNeo2 

JK=14+(0I9-1)72 

S(J)= SY(LOs JK) 

CONTINUE 

DO 85 J=2eNis2 

S(Jd= (8(I41)4S( 5-1) 9/20 

CONT INVE 

DG 87 J=1eNi 

SIT JI=CSY(LKsJ)4+S(05)) 7/20 

CONTINUE 
SJ(N)=0-20 

DO 43 J=1>5N1 

jJJ= 31 

IF(JeEQe-1) GO TO 44 

JJIJ=I—-1 

GO TO 42 

JJJ=JJ 

SYY ={({A1L*SJI( JI) FA2*SY(LK—19JS)4FA3*SY (LK eS JI) FA4*SYIL Ks JIS 
») *AS 

SYN=OME*(SYY-SY(LKeJ))D+SY(LKsJ) 
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DIF=DABS{(SYN—-SY(LK»J))/SSYN 
DI=DMAX1(DI.,DIF) 

SY(LKsJ)=SYN 

CONTINUE 

CALCULATIGN FOR INSULATED SIDE PLATE USING LARGE STEP 
SIZE 

DO 15 [=L0.M 

Tii=i-1 

IFC I2EQ@eM) GO TO 37 

Ti=1+41 

GC TO 38 

Ti=111 

DO 16 J=1laNP 

JJ=J+1 

IF{(Je»EGQ@e21) GO TG 21 

JJIJ=J—-1 

GQ TO 22 

JII=AIS 

IF(12EQ210) GO TO 17 

J1=J 

GO TO 18 

J1i=142*{ J-1) 

SYY={ALL*ESY( ITI, JFA2Z2*SY(T TI s JID + ASBS3*SY(1,j5SS)4A44%SY{ 
xe IsJJJ)) *A55 

GQ 70 19 

SYY=(ALLS*SY( 11,53 JFA22]*SY( ITI s JI) 4¢AS33*SY(15J5JI)4+A44*SY( 
x T,JJJ}) *A55 

SYN=OME* (SYY-SY( I,J) )04S5YT1,/) 
DIF=DABS{(SYN—-SY(1I2J))/SYN 

DI=DMAX1 (DIsDIF) 

SY(IsJ)I=SYN 

CONTINUE 

CONTINUE - 

IF(DI.LE »2ERRe-OR»AKKe»GEsIiT) GO TO 60 
KK=KK+1 

IF(KK.2EQ2200) GO TG 61 

GD TO 70 

DO 63 1I=1sLK 

WRITE(63,67) (SY(14J)0 J=1lseNo2) 
CONT INUVE 

DO 64 I=L0.M 

WRITE16s67) (SY(1s.J) 2J=12NL) 
CONTINUE 

GO TO 70 

WRITE(6262) Kk>5DI 
FORMAT(*°0%5,5Xe"NG GF ITERATION=9 pp I3//5X9"%MAX DIF =' sFi1 
c 0.26) 

WRITE(6352) 
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52 FORMAT(*©*,30X»s *POTENTIAL DISTRIBUTION?) 
DO 65 f=1.4Lk 
WRITE(6,67) (SY(IsJ)sJ=13N22) 
65 CONTINUE 
DO 68 1=i0.M 
WRITE(6367) (SY(1sJ)sJ=15NL) 
67 FORMAT('0*,2X,26F623) 
68 CONTINUE 
N2=N-2 
Ck#¥*x** CALCULATE EZ FRCM SY VALUES USING S-POINT MOLECULE 
DOD 105 I=isLlk 
EZ(%91)9=(-6*SY(155)94+32*SY(154)-—72*SY( 15 3)4+96%*SY(142)-5 
*® O*SY(1,1))7 
14-—24*DyY) 
EZ(%.2)=(SY(1,5)-6*SY(1,4)418*SY(1.3)-10*S5SY(1,2)-3%*SY{ 
aK 131) )/{-12* 
iDdY) 
EZ(CIsNijJ=(SYC1I.N1—3)-6*SY(1I,N1—2 94+1 B*SY(I,Ni-1)—-10*SY{ 
K IsN1)-3*SY(I 
1o»N14+1)97{ 12%*DY) 
EZCIsNI=—-(50*SY(1,N)-96¥*SY(1.N1)4+72 *SY(1,Ni-—1)-32*SY 
x TI,Ni-2)+ 
16*SYCIsN1-3))/( 24*DY) 
DQ 107 J=3»N2 
EZ({isJ)= —-(SY(1Is J—2)-8*SY(1s5—-1)4+8*SY(1,54+1)-SY(IsJ5t+2) 
% I/(12*DY) 
107 CONTINUE 
105 CONTINUE 
WRITE(6.103) 
103 FORMAT(*014,;30Xs *EZ DISTRIBUTIGN® ) 
DO 110 I=1sLKk 
WRITE(639112) (EZ(1sJ) eJ=19N22) 
112 FORMAT(* 04 ,2Xs,21F 8.23) 
110 CONTINUE 
DO 158 fT=1sLKe2 
PUNCH 257s LEZT1sJ)sJ=1+Ns2) 
158 CONTINUE 
NM=NL—-2 
DO 108 I[=LO.M 
EZ(151)=(-6¥*SY(1,5)9432*SY{(154)-72*SY(1.3)496*SY(1,2)-5 
* O*SY(I,1))7 
1(-24*DY1) 
EZ(1s2)=(SY(1,5)-6*S5Y(15,4)418*SY(1,3)-10*SY(1.2)-3*SY{ 
*x Is1))/{(-12% 
1D0Y1) 
EZ( I sNP)I=(SY(I»NF-3)-6*SY(1s+NP—-2)418*SY(IsNP—-1 )-10*SY{ 
Ed IsNP)-—3*SY(I 
1e»NP4+1))7( 12*DY1) 
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EZ(TsNLI= ~(50*SY(IsNL)—S6*SY(IsNP)4+72*SY(1I5NP—1)-32%S 
* Y¥YCESNP-—2)+ 
16*SVYCISNP-3))/( 24*DY1) 
DO 109 J=35NM 
EZ(1isJI= —-(SY( 1, J—2)-B8*SY(155-1)4+8*SY( 1s J4+1)-SYT1I5J5+2) 
* y/(T12*DY1) 
109 CONTINUE 
108 CONTINUE 
DO 115 I=L0.M 
WRITE(63116) (E&Z2(155)5 J=1.NL) 
116 FORMAT(*0'32X%s21F8.3) 
115 CONTINVE 
DO 159 1=LO.™M 
PUNCH 257 s(EZ{L 1s s)sJ=1sNL) 
159 CONTINUE 
Ck#e*e*k CALCULATE EX FROM SY VALUES USING 5S-PGINT MOLECULE 
M2=M—-2 
LT=LK-2 
DD 120 J=14N 
EX(1sJ)I= (-6*SY(55I)432*SY(405)-72*SY(3.I5)4+96*SY( 295 )- 
* SO0*5YT1./) 
1)7{(-24*DX) 
EX(2sJ)I=(SY(5sI)-6¥SYV(455)418*SV(355I)-10*SY(25I)-—3*SY{ 
& 1,J))7{-12% 
1*DX) 
EX{LKKsJ)= (SY(LK—-49J)—-6¥*¥SY{LK—32J)4+18*SY(LK—2s J)-10%5 
* Y{LK-1,J) 
1—-3*SY(LK»J))/(12*DX) 
EX(LK», J)=—-(50*SY(LK+sJ)—-9E6*SY{ILK—-1 J) 4+72*SY(LK—25 J )—32* 
* SY{ILK-32J5)+ 
1 6*SY(LK-45J5))/(24%*DX) 
DO 125 1=3siT 
EX(1,J)=-(SY(1-2545)-8*SY(I-15J5)+8*SY(14+1,/3)-SY(1+2,J)) 
* /(12*DX) 
125 CONTINUE 
120 CONTINUE 
WRITE(62113) 
EX(8s1)=020 
EX(9s1)=0-0 


EX(1121)3= —(SY(9.13)-8*SY{(10,1)48*SY(1221)-SY(13,1))741 
* 2keDX) 
Ex(1001)= —(SY{(8s1)-—8*SY(9.1)4+8*SY(11.1)-SY(1231))/ 12 
x %*DX) 


113 FORMAT(°0%.,30X»e ‘EX DISTRIBUTIN® ) 
DO 130 {=1,LK 
WRITE(69114) (EX(1IsJ)e2J=19N 32) 
114 FORMAT(* 0% s2Xs21F8e23) 
130 CONTINUE 


; ; 


is2 
150 


155 
154 
Ce KKK 


191 


163 


168 
160 


i92 


169 


{CG =LK+1 
L2=LO+1 
L3=L2+1 
0DO 150 J=1sNL 


we i22 


EX(LOsJ)=(-6*SVY(LO+455)432*SY(LO4+35,J5)—-72*SY(LO42 45) 4+96 


* *SY(LO+15/) 
1-S0*SY(LO.J))/(-24*DX1) 


EX(L2sJ)I= (SY(L2435 J) )—-6*SY(L24+ 20J)418*SY(L24+155)-10*SY 


* (L2sJ)-3* 
iSY(L2—-15J4))/(-12*DX1 ) 


EX(M1,J)=(SY(M1 -—35J)-6*SY(M1—255)4+18*SY(MI—1»5J)I-10*SY 


* (M1,J)-3*SY 
1€M5J))/7(012*DX1) 


EX(MsJ)= —(€50*SY(MsJ)—96*SY(M—155)4+72*SY(M —25J5)-32%SY 


* (M—-35J5)+ 
1 6*SY(M—-4,5))/(24*DX1) 
DG 152 I=L35M2 


EX{1,J)=—-(SY(1-22.45)-8¥*SYUI-1,3)48*SY(141.,J5)-SY(1+2sJ)) 


* S(12*DX1) 

CONTINUE 

CONTINUE 

DQ 154 I=LO0..M 

WRITE(64155) (EX{1.j5)sJ=1,NL) 
FORMAT(? 0% 52X%s21F823) 

CONT INUVE 

CALCULATE CURRENT FLCW AT THE CENTER 
WRITE(62191) 

FORMAT(*°O0%s5Xs *CUR IN SMALL STEP SIZE AREA*) 
DO 160 J=i»Ne2 

DO 163 I=1,LK 

CUCTI=EZ( Ts JI) 

CGNT INVE 

CALL DOGSF (DXsCU»ZZsLK) 

WRITE(62168) (ZZ(1),T=1silkKks2) 
FORMAT{L //25F10-6) 

CONT INVE 

ND=17 

WRITE(62192) 

FORMAT(*0%,5Xs *CUR IN LARGE STEP SIZE AREA*) 
JJ=1 

DO 167 J=lsNL 

I=LK 

CUC1)= EZCLKsJJ) 

DD 169 I=LOo™M 

Ji=I1-LO+2 

CU(J1)= EZ(1I+J) 

CONT INUE 

CALL DGQSF (DX1eCUsZZs»ND) 
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167 
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WRITE(6s199) (ZZ(1),1=i,ND) 

FORMAT (//17F1026) 

JJ=JJS+2 

CONT INUVE 

CALCULATE GRADIENT OF JX IN Z-DIRECTION) 

DOG 205 f=1,tk 
CX(1,1)=(-6%EX(1,5)432*EX(154)-72*E X11. 3)94+96*EX(152)-5 
* O*EX(1,1))/7 

if 24*DY) 
CX(1,23=(EX(1.5)-6%*EX(1,4)418*EX(1.3)-10*EX(122)-3*EX( 
a 131) )7{ 12k 

1DY) 
CXCIsN1I=(EX(12N1-3)-6*EX(IsN1I-2)4+18*EX(1sN1—1 )-10*EX{ 
5 IsN1)-3*EX(I 

12N14#1))97(-12*DY) 

CXC(I»Nd= (SO*EXCIsNI-SE*EXLI2N1)472 *EX(I,NI1-—1)-—32*EX 
*x (I.Ni-2)+ 

16*EX(IsN1-3))/7( 24*DY) 

DG 207 J=33;N2 

CX(1i2J)= (EXC 1s J—2 )-B8*EX(15JIJ—1)4+8* EX( 1.341 )-EX( 1.542) 
* 7(12*DY) 

CONTINUE 

CONTINUE 

CX(8s2)= (-6*EX( 696 )+32*EX{( 855)—-72KEX( 894)4+96%EX(823)- 
Es SO*#EX(8s2)) 

1/(24*DY) 

CX(9,2)= (-6¥*EX(926)432*EX( 955 )—-72*EX(934)4+96*EX(95 3)- 
* SO*EX(9s2)) 

l1/(24*xdY) 

WRITE(6s203) 

FORMAT(*04,30X»s "GRAD OF JX IN Z DIRECTION®) 

DO 210 fT=1sLk 

WRITE(63112) (CX( 155) sJ=19No2) 

CONT INUE 

DO 256 I=1.LKo2 

PUNCH 2572(CX(1IsJ)9J=15Ns2) 

FORMAT(3(8F1025/)) 

CONT INVE 

00 208 [=L0.M 
CX(151)9=(-6*EX(1,5)432*EX(154)-72*EX(143)4+96*EX(142)-5 
* O*EX(141))7 

1( 24*DY1) 
CX(C1Ts2)=(EX0125)-6¥*EX(124)4+18*EX( 1, 3)-1LO*EX(1,2)-3*EX( 
* Isl) )/7( 12% 

1D0Y1) 

CXC I sNPI=(EX( I sNP—-3)-6¥%¥EX( 1 »NP-2)4+18*EX(ITsNP-1 )-10*EX( 
% IsNP)-3¥*EX({I 

1 »NP4+1))/{(-12*DY1) 
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124 


CXC IsNLD= (SO*EX( 1,NLI-S96*EX(IsNPI+72*EX( 1 sNP—1)-32*E 
X( I»NP-2)+ 

6*EXTI»sNP-3))/{ 24*DY1) 

DG 209 J=35NM 

CX(IsJ)= (CEXCIs J-2)-8*EX(1,jJ-1)4+8*EX( 15541 )-EX( 12542) 
Y/CL!]*OY1) 


209 CONTINUE 
208 CONTINUE 


* 
1 


x 
1 


* 
1 


* 
1 
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DO 215 I=LO,M 

WRITE(60116) (CX(1IsJ)s J=1sNL) 

CONTINUE 

DO 260 I=LO.M 

PUNCH 258s{(CX(IsJ}sJ=1sNL) 

FORMAT( 3(8F1025/)) 

CONT INUVE 

CALCULATE GRADIENT OF JZ IN X —-DIRECTION) 

DG 220 J=1,WN 

CZ(1sJIJ= (-6%*EZ(59I5)432*EZ1 43/5 )—72*EZ(355)+96%EZ1235)- 
5S0*EZ(1.J5) 

IS { 24%*DX) 

CZ(20I)=1(EZ(55I5)-6¥EZ145 JIF1B*XEZ( 35 J)—-10*EZ( 25 J )—-3*EZ{( 
l»J)})/( 12% 

*DX) 

CZ(LKKs JJ=—-(EZ(LK—-4,.4 J)-6¥*EZ(LK—-33s45) +1 8*EZ(LK—2s J)-10*E 
Z{LK—-13J) 

~3*EZ(LKeJd))/(0124DX) 

CZIL Ks J) = (SO*XEZ{IL Ks J)—-GO*XEZ(LK—-15J )4+72*EZ{LK—22 J )-—32* 
EZ(LK—3.J)+ 
6* EZ(ILK—-4545))/(24*DX) 

DG 225 [=3sLT 

CZ(12eJ)= (EZ 1-2.J5)-8FEZTI-12J5) +8*EZ(141,5)-EZ(14+22J)) 
4CT12*DX) 


225 CONTINUE 


220 


Pot Be 


230 


270 


* 
1 


x 


CONTINUE 

WRITE(6.213) 

FORMAT(40%s.30Xs *GRAD GCF JZ IN X DIRECTION® ) 

DO 230 I=15,LK 

WRITE(62114) (CZ(1Tsd)a J=1sNo2) 

CONTINUE 

DO 270 I=1»sLKs2 

PUNCH 2579 (CZ(1sJ)sJ=1sNs2) 

CONTINUE 

DO 250 J=1sNL 

CZ(LOs J) =( -6¥*EZ(LO4+45I5)432*EZ(LO4+3, J)-72*EZ(ILGO+2,5)+96 
*XEZ(LO+155) 

—-SO*EZ(LCsJ)II/( 24*DX1) 

CZ{(L2.J5)= (EZ(L24355)-6¥*EZ(L24255)418*EZ(L24+1,.5)-10*EZ 
(L2sJ)—-3* 
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LEZ(L2—-1245))7( 12*DX1) 
CZ(M15JI=(EZ(IM1 —3sJ)-6¥EZ(M1-2sJ)+18EZ(M1—-1,/5)-10¥EZ 
* (M1,J)—-3*SY 
1{MsJ))/(—-12*DX1) 
CZ(MeJ)= (50* EZ (Ms J )—-964EZ(M—155)4+72*EZ(M —22J5)-32*EZ 
* C(M—353J5)+ 
1 6*EZ(M—-445))/(24*DX1) 
DO 252 I1=L3.M2 
CZ(1sJ)= (EZ(1I-2.J5)-8*EZ(1-15J5)4+8*EZT1I4+1,j))-EZ(14+23J)) 
* S/(12*DXK1) 
252 CGNTINUE 
250 CONTINUE 
DG 254 I1=L0,M 
WRITE(63155) (CZ(1sJ)sJ=isNL) 
254 CONTINUE 
DO 265 I=LOomM 
PUNCH 258s {(CZ( IJ) sJ=1+sNL) 
265 CONTINUE 


STOP 
END 
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* IN TRE SIXTH COLUMN INDICATES 
CONTINUATICN FROM PREVICUS LINE 
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FORTRAN PROGRAM FOR CHAPTER 5 

THIS PROGRAME CALCULATE PRESSURE DOISTRIBUTIGN FOR 

FINITE WIOTH PARALLEL PLATE SLIDER BEARING 

THE PRESSURE EQUATICN IN FINITE DIFFERENCE FORM IS 

SOLVED USING POINT SUCCESSIVE OVERRELAXTION METHOD 

SUBROUTINES FROM IBM-SYSTEM/360 SCIENTIFIC SUBROUTINE 

PACKAGE ARE USED IN THIS PROGRAM 

DOUBLE PRECISION A0(41541) 2sA1sA2sA3sA4 a AS sHART(5)2DXsD 
ZsaSLOP>, 


1P{41241)2PP.,PN sOMEZERRsDIsDIF.Z(41) 52271 41)5Z7Z1f41).TCX(4 


1s41)s 


2027041541) »EZ(141 141) 5QX1(41)5QX2141) sPGR141241) 


DATA FOR CALCULATION 
M=41 

N=41]1 

NN=21 

M1i=M—1 

Ni=N-1 

DX=02025 

DZ=0.2025 

OME=1.4.5 

ERR=0.00005 

IT=300 

STARTING VALUES GCF PRESSURE AT GRID POINTS 
DO 22 I=15M 

DO 32 J=1.NN 
P(Iy9J)=020 

CONTINUE 

CONTINUE 

READ HARTMANN NUMBER 
DO 15 L=1,51 
READ(Se11) HART(L) 
FORMAT( F1026) 
WRITE(6912) HART{L) 
FORMAT(*0%+s5Xs "HARTMANN NUMBER="'4.F1026) 
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110 
10 


20 


118 
120 


122 
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27 
25 
CHK 


45 


47 


49 
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READ CURRENT CENSITY GRADIENTS AS QBTAINED FROM THE 

SOLUTIGN OF POTENTIAL EGUAT IGN 

DO 10 {[=1,41 

READ(53110) {CXC1+sJ)sJ=1221) 

FORMAT{®3¢ 8F1025/))} 

CONT INUVE 

DO 20 f=1,41 

READ(52110) {(CZ{1IsJ)+J=1321) 

CONTINUE 

DO 120 1=1,41 

WRITE(6,118) (CX(I,JI),J3=1221) 

FORMAT(2X/21F1i0.25) 

CONT INUVE 

DO i22 {=1341 

WRITE(Gs118) (CZ(IsJ)5J=1,21) 

CONT INVE 

CALCULATE COEFFICIENTS OF FINITE DIFFERENCEFORM EGQe 

A1=1./{(DX**2) 

A2=12e/fDX**2) 

AB=10e/(DZ**2) 

A4=16/1(0Z**2) 

AS5S= (0D0X4#DZ )**2)/1724(DX**24DZ* 2) ) 

DO 25 I1=25M1 

DO 27 J=2aNN 

AG (@19J)=-HART(LI*(CX(1isI5)4CZ1155)) 

CONTINUE 

CONT ENUE 

START ITERATING 

KK=0 

DI=0.20 

DO 35 iI=2sM1 

{i=T +1 

Tii-iI-1i 

DO 40 J=23NN 

JJ=J3+1 

JJ jJ=JI~-1 

IF(JeEGe21) GO TO 47 

PP=( AOL ITsJ) FAL FP CII, J)FA2*PLILIs J) 4A3*PL IIS) FA4*P( I 25 
JJIJI#XAS 

GO TO 49 

PP={ AOL IT »sJIFAL¥PC IIs J) +A2*PL ILI s J+ A3¥P (ITs J ISJD+A4#P TI 
JJJ)) *AS 

PN=DME*(PP-—-P{(1I,J)) +P(1IsJ) 

DIF= DABS{PN—-P(I5J))/PN 

DI=OMAX! (DI,DIF) 

P(1IsJ)=PN 

CONT INVE 

IF (DIeLEeERReOReKKeGE2IT) GO TO 50 
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KK>KK4+1 
GO TG 45 
50 WRITE 6255) 
55 FORMAT(4*0?#,25Xs *PRESSURE DISTRIBUTION® } 
WRITE(6256) KKkK,sDI 
56 FORMAT("0*.5Xs *NO OF ITERATIGN=* s137/95X%s*MAX DIFF*,F1 
* 0426) 
DO 60 J=1.NN 
WRITE(6,58) (PC1sJ})sI=15Ms2) 
58 FORMAT(*0*.2Xs21F 1023) 
60 CONTINUE 
CeKeEX CALCULATE LOAD CAPACITY 
ND=41 
DD 70 J=2sNN 
00 75 1=1,M 
Z{UIP=P(1T.J5) 
75 CGONTINVE 
CALL DOSFL{DXKeZsZZaND) 
WRITE(6.80) €ZZUT1),T=1eMs2) 
BO FORMAT(I//2Xa21F623) 
Zi (J3=ZZUND) 
70 CONTINUE 
Z1(439=0.0 
CALL OGSF {(DZsZ1eZZ2NN) 
WRITE(6:85) ZZ(NN} 
BS FORMAT(S/5Xse *LOAD="sF1046) 
15 CONTINUE 
STOP 
END 
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